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ABSTRACT 


The development of systems such as acoustic levitation chambers will allow for the 
positioning and manipulation of material samples (drops) in a microgravity environment. This 
provides the capability for fundamental studies in droplet dynamics as well as containerless 
processing work. Such systems utilize acoustic radiation pressure forces to position or to further 
manipulate (e.g., oscillate) the sample. The primary objective of this report was to determine 
the effect of a viscous acoustic field/tangential radiation pressure forcing on drop oscillations. 
To this end, the viscous acoustic field is determined. Modified (forced) hydrodynamic field 
equations which result from a consistent perturbation expansion scheme are solved. This is done 
in the separate cases of an unmodulated and a modulated acoustic field. The effect of the 
tangential radiation stress on the hydrodynamic field (drop oscillations) is found to manifest as 
a correction to the velocity field in a sublayer region near the drop/host interface. Moreover, 
the forcing due to the radiation pressure vector at the interface is modified by inclusion of 
tangential stresses. 
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I. INTRODUCTION 


IA. BACKGROUND 

The development of acoustic levitation systems has provided a technology which will 
enable containerless processing in a microgravity environment. Such systems also provide the 
capability to undertake fundamental studies of droplet dynamics. Acoustic levitation devices 
utilize radiation pressure forces to position the sample away from the chamber walls. Specific 
types of systems are the single axis system (Barmatz, 1981) as well as a three-axes design (Wang 
et. al., 1984). In the case of a three-axes system, acoustic drivers (speakers) centered in the 
three orthogonal sides of a parallel piped chamber are driven at its resonant frequency. Liquid 
drops can then be positioned in the region in which the pressure is at a minimum (i.e., the wave 
pressure nodes). The acoustic levitation system provides for more capabilities than merely that 
of drop positioning. In addition, the drop can be made to rotate via phase lagging two of the 
acoustic waves. Also, drop oscillations can be induced via frequency modulation of an acoustic 
wave. 

The quadrupole resonance of simple drops was investigated by Marston and Apfel (1980) 
in an experimental study. Modulated acoustic radiation pressure provided the driving force. 
The experimental work of Trinh, Zwem, and Wang (1982) studied the small amplitude 
oscillation and decay of a free (non-driven) drop. Furthermore, large-amplitude drop shape 
oscillations have been experimentally investigated for the free and forced cases (Trinh et. al., 
1982). Both drop oscillations and break-up were studied in a visualization experiment by 
Marston and Goosby (1985). All these experiments were performed in a lg gravity field. The 
drop itself was surrounded by an immiscible host fluid, with positioning accomplished using 
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acoustic radiation forces. 


Analytical and numerical investigations have focused primarily on free drop oscillations 
and decay. Miller and Scriven (1968) addressed the problem of normal mode oscillations of a 
viscous liquid drop immersed in a viscous host medium, and derived analytical expressions for 
limit cases. More recently, Prosperetti (1980) completed a numerical study of free viscous drop 
oscillations and decay for a range of values of host medium viscosities. Both these 
investigations involved a linear analysis. Also, there has been some recent work concerning the 
nonlinear oscillation of inviscid drops (Tsamapoulos & Brown, 1983, Natarajan & Brown, 
1987). 

The work of Marston (1980) was an analytical investigation into the acoustically forced 
fluid drop problem. It is the projection of the radiation stress tensor onto the surface of a drop 
in order to form the radiation pressure vector which actually accounts for the forcing terms. 
Marston calculated the radial component of the radiation pressure vector for a specific limit case; 
that in which the acoustic field is taken to be irrotational and the correct tangential acoustic 
boundary conditions at the drop/host medium interface were not enforced. Thus, no tangential 
radiation stresses could contribute. The drop was forced at the boundary/interface only. There 
was no forcing of the Navier-Stokes equations which govern the behavior of the hydrodynamic 
field which itself results from the action of the acoustic forcing. The only coupling between the 
hydrodynamics and acoustic fields is through the boundary/interface conditions. 

The incorporation of viscous effects into the acoustic field has ramifications for the full 
hydrodynamic problem in two ways. First, the modification of the acoustic field to include 
viscous effects will allow for the enforcement of the tangential boundary condition on the 
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acoustic field. This will result in the existence of tangential radiation pressure forces. These 
forces then enter into the boundary/interface conditions on the hydrodynamic field. Secondly, 
the full hydrodynamic equations themselves are forced by nonlinear terms relating to the acoustic 
field. The solution of the hydrodynamic field itself is then altered from that given by, say, 
Miller and Scriven (1968). 
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IB. OBJECTIVES 

The primary goal of this work is to determine the effect that incorporation of the 
tangential radiation pressure forces has on the drop oscillations. 

Additional related considerations include: (a) possible effects of a non-axisymmetric 
acoustic field, (b) efficient evaluation of the tangential radiation pressure vector, (c) the role of 
bulk viscosity in the tangential radiation stresses, and (d) extensions to compound drop forcing. 

In order to meet the primary objective, it is necessary to include the effects of viscosity 
in the acoustic field representation. This is done in Section III of the report. Since the viscous 
acoustic field will affect the hydrodynamic field both through the boundary/interface conditions 
and via a modification of the hydrodynamic field equations themselves, these topics must be 
addressed. The determination of the tangential radiation pressure vector which contributes to 
the boundary/interface conditions is done in Section IV. The modified hydrodynamic field is 
investigated in Section V for the case of forcing by an unmodulated acoustic standing wave field 
and in Section VI for the case in which the acoustic standing wave field is modulated. 

The secondary objectives are addressed in the appropriate sections (see Table of 
Contents). Moreover, supporting material in the Appendices serves to further elucidate these 
topics. 

Conclusions are presented in Section VII. 
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IL ACOUSTIC FORCING MECHANISM AND GOVERNING EQUATIONS 


IIA. ACOUSTIC FORCING MECHANISM 

Certain definitions and an associated discussion will be useful, as they pertain to concepts 
and quantities which will be calculated/manipulated in subsequent sections. 

The acoustic radiation stress tensor is given by 


MaJ “ P (HA.l) 

with 


P = i<P‘> -(■£-)< V A > 


(HA. 2) 


Adiabatic compressibility is given by P . The brackets indicate averaging over an acoustic 
period. The projection of the radiation stress onto the drop surface is given by 



( / C - TC ) 

+ (TC - TC) + (T/r - T£)e f (iia.3) 


Superscripts indicate the regions interior and exterior to the drop’s surface. The tangential 
component of the radiation pressure vector is given by the and e* components of equation 
(IIA. 3). Note the velocity and pressure fields indicate acoustic quantities. 

The drop oscillation is induced through frequency modulation of the acoustic wave. 
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Given an acoustic frequency, , a second wave with frequency u 2 is required, such that 


the drop frequency of oscillation is given by the difference of those of the acoustic waves. That 
is, ^ ( drop ) ~ u r 03 2 • The carrier frequency is given by the average of the two acoustic 

frequencies. Note that the wavelength of the acoustic carrier wave should be such that it couples 
reasonably to the drop diameter. The representation of the acoustic radiation pressure vector 
is best done through the use of spherical harmonics. That is, 



P r 


The projection is 


FW + prwstui*-^ 


(IIA.4a) 


(pr)^ -v (pr)^ 


ST MIC 


Jl 




1 


(HA.4b) 


and similarly 


, xTAMG, 
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(IIA.4C) 


Note that <2Q-sin0d0c?4> . The static part of the radiation stress vector as well as the 
oscillating contribution (which oscillates at the drop frequency) are given. If the acoustic wave 
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is not modulated, there will NOT be an oscillating forcing term. Note that the tangential 
radiation pressure appears in an odd operator form. This is due to the fact that the unforced 
drop oscillations were originally determined (Miller and Scriven, 1968) using this form. 

For further discussion of the radiation stress as related to drops, see Marston et. al., 

(1982). 
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IIB. CASE OF NON-AXISYMMETRIC FORGING 


This section pertains to a qualitative investigation into the possible forcing of shear waves 
(in the hydrodynamic field) due to the acoustic field. As such, the acoustic field will be 
assumed to be viscous, and to generate non-zero tangential radiation stresses at the interface. 
As effects can be seen even if the coupling between the acoustic and hydrodynamic fields is 
restricted to the boundary conditions, this will be the taken to be the case. That is, no forcing 
of the governing equations themselves will be taken into account. Again, this is making the 
inviscid acoustic approximation in the governing equations of the hydrodynamic field, but not 
in the boundary conditions. This introduces the simplification that the general solutions to the 
unforced hydrodynamic equations are those given be Miller and Scriven (1968). 

The task of this section is to investigate the result of boundary/interface forcing on the 
behavior of these solutions. The boundary conditions which must hold at the drop/host interface 
are (1) the kinematic condition, (2) the stresses are balanced, and (3) the velocity components 
are continuous. Both the second and third conditions generate three equations. Moreover, all 
the boundary conditions can be expressed using only the radial component of the velocity and 
vorticity fields. This can be done via operations of taking the surface divergence of the 
tangential balances in conditions (2) and (3), and of taking the radial component of the surface 
curl of these same conditions. A spherical coordinate system is used. 

For clarity, the radial velocity (fl r ) and vorticity (6 r ) components are listed below 




+ 





(HB.la) 
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(KB. lb) 




The functions are spherical Bessel and Hankel functions of the first kind. The caret indicates 
the hydrodynamic velocity field. The ( 0 , <|> ) are spherical harmonics . 

Note that the "a" and "b" coefficient are determined via application of the boundary 

conditions. In this case, the boundary conditions will involve forcing due to the radiation 

stresses at the interface. The forcing of the shear waves is constructed via performing the 

aforementioned operations involving the surface curl. This yields (at r=R) 


b < (-£ji(U)) - k - 0 


(KB.2a) 
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(UB.2b) 


The exponential term indicates the phase. Without the forcing term on the right hand side of 

equation (IIB.2b), this system will reduce to two homogeneous equations in two unknowns. The 

resulting dispersion relation describes the behavior of shear waves and is well known. 

Physically, it describes waves which are purely rotational and decaying. The presence of the 

forcing term is due itself to the components of the tangential radiation stress vector. Note that 
if axisymmetry is assumed, and pr 4 is zero, then there is no forcing. However, if the 
axisymmetry of an acoustic field is not maintained in, say, an experiment, then the forcing of 

hydrodynamic shear waves appears possible. Of course, a full investigation would require that 

the hydrodynamic field equations themselves be forced in addition to the boundary/interface 

forcing described above. 
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EC. GOVERNING EQUATIONS 

The governing equations of those of fluid mechanics; the conservation of mass, 
momentum, and energy equations. Since the fluid (acoustic field and hydrodynamic field) is 
considered to be isothermal, the energy equation will simplify considerably, and result in a 
relationship between the pressure and density fields. 

The general form of the governing equations will include contributions due to 
compressibility. Although the hydrodynamic field is considered incompressible, the acoustic 
field clearly cannot be. 

The ambient mean pressure and density are constant. The acoustic field (or sound field) 
can be viewed as a perturbation of the ambient environment. A consistent perturbation 
expansion scheme will be developed whereby the hydrodynamic field will arise naturally (at 
second order), forced by terms which are quadratic in the acoustic field quantities. 

The governing equations must be satisfied in the regions interior and exterior to the drop. 
In dimensional form they are given by: 

(Conservation of Mass) 

r\ - 

— + k. • **■ P^’U; - O 

(HC.l) 

(Conservation of Momentum) 

f(-l- 

= -Vp + + (^-+yU ButK )V(V- U,) (DC2) 
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(Conservation of Energy/Equation of State) 


|p = 1 p (nc.3) 

The variables p, u , and p denote the pressure, velocity and density fields. The material 

properties of dynamic viscosity, p , and bulk viscosity, V- BULX , are constants. Of course, the 
values will be different according as to which fluids form the drop and host media. The 

conservation of momentum equation is the compressible Navier - Stokes equation in which no 
body forces are present. In equation (IIC.3), c Q refers to the speed of sound. 
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m. VISCOUS ACOUSTIC FIELD 


m A. DISCUSSION OF APPROACH: VELOCITY DECOMPOSITION METHOD VS. 
METHOD OF COMPOSITE EXPANSIONS 

Attention is restricted to the case in which the acoustic field is axi symmetric. In the case 

of the velocity decomposition approach, the acoustic velocity field is comprised of a solenoidal 

and an irrotational contribution. That is, 

V = V X A + V<J> (ni.1) 

Moreover, A is of the form (0,0, ijr) , with i|r a scalar function. Equations can be developed 

for <J) and ijr . The details of such a development, and the resulting solutions, are found in 
Appendix I. That is, for the axisymmetric case, the problem of determining the acoustic field 

which incorporates the effects of viscosity is analytically very tractable. 

In particular, note that arguments (of spherical Bessel functions) arise in the form of 
^SQRlJito/v^jR , with o the acoustic frequency and R the radius of the drop. Specifically, 

these arise in the solution for ijr . (Of course, ijr itself contributes to the boundary conditions 
at the drop/host interface.) The acoustic frequency is on the order of hundreds of kilohertz, and 

the drop radius is, in general, several millimeters. Moreover, it is necessary to allow for 
reasonable values of \ g , so it is doubtful that realistic values of the kinematic viscosity would 
serve to decrease the magnitude of the argument. 

This leads to difficulties in the numerical evaluation of the unknown coefficients (listed 

in the boundary/interface conditions given by equations [Al.lla - AI. 1 Id]). 

An alternate approach which avoids these difficulties is to utilize a boundary layer 

formulation. Moreover, this approach has the capability to elucidate the flow field structure in 

the region near the host/drop interface. 
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It is this second approach which will be employed. However, note that the first approach 
has been developed, also, and is discussed in detail in Appendix I. 

The boundary layer approach which is to be used is termed the method of composite 
expansions (see Nayfeh, 1973). The details of this method as they pertain to this problem are 
found in succeeding sections. 
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Ill B. SOLUTION OF THE VISCOUS ACOUSTIC FIELD VIA THE METHOD OF 
COMPOSITE EXPANSIONS 

The governing equations are those listed in Section II C. They will be 
nondimensionalized with respect to acoustic field variables. Let 

(“’* x = )<. 

(<£) v=v ; p‘f ~ = p (nni.i) 

The velocity scale is given by the sound speed in the outer medium. The acoustic frequency is 
used to define the time scale. The length scale follows from these two quantities. Density is 
nondimensionalized with respect to the ambient density in the outer medium. Tildes indicate 
nondimensional quantities. The nondimensionalization scheme is the same in both regions, 
interior and exterior to the drop. 

The method of composite expansions is a generalization of the method of matched 
asymptotic expansions which has been used successfully in boundary layer applications (Nayfeh, 
1973). The basis idea is that the viscous effects on the acoustic field are of primary importance 
in a region near the drop/host interface. The scale of this region is proportional to SQRT (V/cu), 
with v the kinematic viscosity and w the frequency of the acoustic wave. In levitation systems, 
this quantity is quite small due to the relatively large magnitude of the acoustic frequency. The 
choice of the acoustic frequency is driven by the selection of the acoustic wavelength to be on 
the order of the drop dimension. 

The full nonlinear nondimensional governing equations of compressible isothermal fluid 
dynamics are to be expanded in a small parameter, 6 . In this analysis, 5 is a formal 
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parameter. It is interpreted as the ratio of the natural frequency of oscillation of the drop to that 

of the acoustic frequency. Another small parameter arises naturally in this ’’boundary/interface 

layer" approach. This second small paramter, e , is equal to SQRT (1 /Re) , with Re a 
Reynolds type number which is defined later in this section. 

It is possible to relate the formal structure of this analysis to work done by Riley (1967) 

in his consideration of the flow induced by a solid body that is in oscillatory motion in an 

infinite viscous fluid which is otherwise motionless. Of course, the problem being considered 

here is very different. The fluid drop has no translatory motion. Oscillations of the drop would 

be due strictly to the action of a modulated acoustic field. However, a formal comparison may 
be made with Riley’s work. He defined three parameters, R, Rs, and M 2 . These are equivalent 

to (6/e 2 ), (6/e) 2 , and (l/e 2 ) , respectively. The region of interest in the acoustically forced 

drop problem would correspond to the case of M > 1, R, Rs > I . 

Consider the fluid region exterior to the drop (host region). The system of governing 

equations (Section IIC) is nondimensionalized and then is linearized after a perturbation 

expansion in a small parameter, 6 . The base state about which the perturbation occurs is that 
of zero mean motion, and constant ambient pressure and density. The lowest order perturbation 

is the sound field. The tildes are dropped for convenience. The resulting system (to lowest 
order in 6 ) is 


2Pl + V- V,° =0 (HB.2a) 
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(HB.2b) 
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0 s. 
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(mB.2c) 


The superscript "o" denotes the fluid region exterior to the drop. A subscript of "1" indicates 
sound field quantities. At this stage, higher order terms have been neglected. The 
nondimensional parameters which appear are 


t _ 


T* - Km*/)C 


Rft (cf/ldfit) 

In the fluid region interior to the drop, the governing equations (to lowest order in 
small expansion parameter 8 ) are 


Ifl + ft v • v " = o 

it r 


(mB.3a) 


av, 4 ' = 


ht 


(f)v + 


+ Ik) ) (mB - 3b) 
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(nm.3c) 




with a - p - p[/po , and x‘ - p^u/pjj • The system of equations (IIIB.3a-c) in the 

region interior to the drop (and (IIIB.2a-c) exterior to the drop) are subjected to a "boundary 

layer" type analysis. That is, in the interior of the drop, the fluid will be considered to be 

comprised of an "inner" layer and an "outer" region. Note both the "inner" and "outer" regions 

lie within the drop. In the "inner" region, the effects of viscosity are important, and in the 

"outer" region they are absent. This analysis itself is termed the method of composite 

expansions. The same methodology is employed in the (host) fluid exterior to the drop. See 

Figure IHB.l. 

The subscript "1" which heretofore has referred to sound field quantities is now dropped. 
Subscripts will now refer to inner and outer expansion quantities . 


Application of Method of Composite Expansions to Linearized System of Equations 

In the method of composite expansions, the dependent variables are expressed in terms 
of functions having different independent variables (in the e r direction). One is a stretched 
variable. In the outer region, the variable is given, of course, by r. In the inner region, the 

independent variable is given by 5 (if interior to the fluid drop) or C (if exterior to the fluid 
drop). The small parameter in the expansion which is part of the method of composite 

expansions is e . This small parameter is taken to be larger than the small 

parameter 6 (found in the 6 expansion). That is, this method is valid for the case in which 
the nonlinear terms arising at next order are smaller than the viscous contribution. 
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Although this approach is valid when the acoustic field is non-axisymmetric, the equation 
development will be done for the less algebraically involved case (i.e., axisymmetric). 

Consider the fluid region interior to the drop. Let r-R-e £. The constant 

R - au/c° , and "a" is the nondimensional drop radius. The small parameter, e , is taken 
to be sqrt(l/Re). Rewriting system (IIIB.3a-c) in terms of the inner layer variable, £ , yields 
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The dependent variable q‘ E {v r ‘, v„', p‘, p'} is expanded in the inner layer region as 
follows: 

V = ( + V(w>) 

+ £ ( $m* - 5 %»1s + 

♦ H.*T. (,Iffi - 5) 

The notation nA "' indicates the outer region dependent variable. It is expanded in a Taylor series 

expansion in terms of the inner region variable. The ” ~ " denotes the inner region dependent 

variable. Superscripts "i" indicate that these quantities are defined in the fluid region interior 
to the drop. The inner region functions must decay to zero as 5 - 00 • 
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The functions which are dependent upon the outer region variable, r, can be obtained via 

solution of outer region system. To lowest order in e , this system is simply that which 

represents the inviscid acoustic standing wave field. For completeness, these equations and the 
resulting 

solution are listed below 


+ p.V 1 sO 

n 


(11 IB. 6a) 


■jT 



(IIIB.6b) 




A solution is obtained in <f>‘ , which is 
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j, is a spherical Bessel function. The coefficient is 
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with h the distance from the acoustic velocity nodal plane to the center of the drop. If this 

distance is zero, then t , is zero for odd 1. The coefficient t , should not be confused with 
the small expansion parameter 5. P, is a Legendre polynomial. The other dependent variables 

in the outer region (still interior to the fluid drop) can be constructed from knowledge of <$>' . 

Note that the unknown coefficient (at each value of 1) is given by aj . 

Substitution of the expansions (IIIB.5) into the system (IIIA.4a-c) yields 
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Note that 



n 



f K 


from the outer region analysis (that is, from the solution of the inviscid acoustic field interior 

to the drop, which is then expanded in a Taylor series expansion in terms of £ , the inner 
region variable . This equation is embedded in (IIIB.7c). Similarly, 



e f ( Ktt ) + «• P-iUiK ♦ £ P 1 


0 


which is embedded in (IIIB.7a). Also, 



which is found embedded in (IHB.7b). From the remaining terms in (IUB.7a-c), it can be seen 

that vj| = 0(1) =>^ » = 0(e) => , and p*, p'~ (K® 2 ) -* py Pi • Then, to lowest 

order , the inner region equations (interior to the drop 1 ) are 
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The solution for v/ can be readily obtained. Following that, it is straightforward to 

°o 


obtain vf t , and then p[ . The solution for ^ is given as 
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The determination of the acoustic field quantities in the fluid region exterior to the drop 
(ie., host region) may be done in a similar fashion. An inner variable, C is defined. Let 


r-R+e f. Again, the small parameter e equals sqrt (1/Re). Also, R - aoi/c° , with "a" 

the dimensional radius of the drop. As was the case in the fluid region interior to the drop, the 


dependent variables are to be expanded as 
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(DIB. 10) 


where q° E {v r °, v„°, p°, p°}. The superscript "o" denotes the fluid region exterior to the drop. 


The outer region dependent variables are indicated by " A ", and the inner region dependent 

variables are known by the " Functions in f must decay to zero as ( - °° . Note the 
Taylor series expansion form of the outer region dependent variables are expressed in terms of 

C , the inner variable. 


The functions which are dependent upon (r,0,f) can be obtained by solving the outer 

region system. To lowest order in e , this is the inviscid (scattered) acoustic standing wave 
field. For completeness, the equations and resulting solution are listed 
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with v 0 = V <j£ T> A solution is obtained for <j£j, which is 
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(mB.lld) 


h, (1) is a spherical Bessel function. Both P, and o x are as given previously in this section. The 


unknown coefficient (for each 1) is given by a.° sl . The other dependent variables (in the outer 

region and exterior to the drop) may be generated from knowledge of 4»^ T - 

The system of governing equations exterior to the drop (IHB.2a-c) rewritten in f (inner 
variable), yields, to lowest order in e , 
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The solution for can be readily obtained. Following that, it is straightforward to obtain 

CO ~ o oi O 

v rr and the p 2 . The solution of v 4 is 
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(DIB. 13) 


The solutions for v p , , v fl , p 2 , and p 2 are listed in the following. For 

completeness, ^ and v£ g are included. Note that to lowest order, only the e e component 
of velocity has a contribution in the inner variable (f for f for \g ). It is at higher orders 

that the other components will contribute terms which are strictly functions of the inner variable. 
In order to determine the total contribution at order e , it would be necessary to find ^ and 

v^ . In a similar fashion, in order to determine the total contribution at order e 2 , it would 
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, at© * 0 

be necessary to calculate v r* > v fj * ^ ^ • Moreover, the "outer" region solution must be 

modified at this order. 

These complete sets of higher order contributions are not needed. However, if there 
were to be interest in, say, the function v^ (at order e ), it is noted that simply solving for 

them via the order e equations which derive from (IIB.7a-7c) would not yield a uniformly valid 
expression for . That is, secular terms arise. These can be eliminated via introduction of, 

say, a slow time scale. Such calculations are not required in the execution of this project, and 

so are not presented here. 

It is found that 
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It is stressed that, to 
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To reiterate, the " A " indicate the solutions in the outer regions, which correspond to inviscid 
acoustic field solutions. 
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mC. BOUNDARY/INTERFACE CONDITIONS 


Solutions have been obtained for the inner and outer regions, both interior and exterior 
to the drop. All solutions are finite in their respective regions of validity. The unknowns (for 
each value of 1) are a\, ct° h . 

There are four boundary/interface conditions which must be imposed. Use of these 
conditions will generate a set of algebraic equations (at each 1) which can be solved for the 
unknowns. Since the equations are nonhomogeneous, that is, they are forced by (functions of) 
the incident acoustic wave, this is not an eigenvalue problem. Once the coefficients are known 

the viscous acoustic field, both interior and exterior to the drop, is known. 

The four boundary/interface conditions at r-R 
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The boundary/interface conditions must be applied at r-R , which corresponds to £ = 0 and 
f = 0. Note that in this context r refers to the stresses. Physically, these boundary conditions 

represent the velocity and stress balance in the normal and tangential directions to the interface. 

If the solution of interest had been restricted to just the inviscid acoustic field, only 

conditions (IIIC. la-b) would be applicable. Moreover, r„ would be replaced by pressure. 

To lowest order in e the conditions are 
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In (HIC.2a-b),JE, > 1. If the acoustic field were inviscid only, then the acoustic field variables 

would be known once a\ and a° sl were known (at each 1 value). That is, a system of two 
algebraic nonhomogeneous equations in two unknowns would have to be solved. The forcing 

term is due to the presence of the incident acoustic standing wave field, without which there is 

no levitation system or flow field. 

In the first boundary /interface condition, it is the radial components of the acoustic 
velocity which must balance at the interface (r-R, C-5-0) . Note that this condition involves 

only the contributions from the outer regions (i.e., inviscid solutions), as the viscous correction 
occurs at a higher order. This is also the case for the second boundary/interface condition in 
which the viscous contribution to the normal stresses occur at higher order. 

It is in the tangential acoustic velocity balance at the interface that contributions from the 
"inner" region enter, along with those which represent the inviscid solution. Finally, in the 
tangential stress balance, only the inner region solutions contribute. Therefore, it is tangential 
velocity balance which provides the coupling between the inviscid "outer" system and the 
"inner" region corrections. 

Since viscous effects are to be taken into account, conditions (IIIC.2c-d) also contribute; 
resulting in four equations in four unknowns which must be determined. 

This solution is done numerically using parameters of interest. Such quantities include 
those of viscosity and density inside and outside the drop, as well as the drop radius. Note that 
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the bulk viscosity does not appear (at lowest order) in either the outer or inner region equations. 

Knowledge of the viscous acoustic field is necessary in order to determine the tangential 
stress component of the radiation pressure vector. This is done in the next section. 

Moreover, if the acoustic field is viscous, the acoustic field will couple to the second 
order hydrodynamic field not only through the boundary/interface conditions, but also act as a 
forcing term in the hydrodynamic equations themselves. This is presented in upcoming sections 
for both the unmodulated acoustic standing wave field case and the modulated acoustic standing 
wave field case. 
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HID. COMMENTS ON THE ROLE OF BULK VISCOSITY 


In the previous section, it was shown that part of the contribution to the attenuation of 
the acoustic field is due to the presence of bulk (expansive) viscosity which enters into the 
governing equations. The equations governing the acoustic pressure and velocity fields are 
solved retaining this physical contribution, in general. 

However, in the determination of the viscous correction to the acoustic field via the 
method of composite expansions, it is found that those terms involving the bulk viscosity only 
contribute at higher orders. Thus, in this approach, the bulk viscosity is not relevant. 

If the vector decomposition approach had been utilized, the bulk viscosity would enter 
into the analytical solution. This can be seen in Appendix I (see for example, AI. 6 B). It acts 
to attenuate the acoustic wave. 

However, in order to determine actual numerical values, it is necessary to know physical 
property data. This necessity is not restricted to the value of the bulk viscosity, but also applies 
to ambient medium density, shear viscosity, speed of sound in the medium, etc. 

Of course, for a monatomic gas, the bulk viscosity is zero. Among the empirical 
formulas which give values for n is the classical work of Greenspan (1959) concerning 
attenuation in diatomic gases. He found (using a polynomial fit to the data) 

JU. = /(.(lY* - 1.3** -(Jr-i)/?r ) 

where n (the shear viscosity) is taken to be constant, a is the empirically determined constant, 
7 the ratio of specific heats, and Pr the Prandtl number. For air, a = 0.8903, and for nitrogen, 
a = 0.9104. 
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When a specific physical system is under consideration, it is necessary to have at hand 
the relevant physical property data in order to determine numerical values. However, the effect 
of the physical phenomenon of bulk viscosity (among other contributing factors) is retained in 
the analytical formulation which involves the exact solution (rather than the approximate 
expansion method). 
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rv. TANGENTIAL RADIATION PRESSURE VECTOR 


The calculation of the radial component of the radiation pressure vector was done by 
Marston (1980) for the case in which the acoustic field was considered to be strictly inviscid. 
This topic will not be readdressed in the report. That is because the additional complexity due 
to the viscous contribution will be evident from the full discussion of the tangential contribution 
to the radiation pressure vector. Also, the recalculated radial radiation pressure would be, in 
some sense, a correction to that done for the inviscid case. Moreover, in the case in which the 
acoustic field is considered to be fully inviscid, as was done by Marston, there is no contribution 
whatsoever to the tangential radiation stresses. Thus, no counterpart to the following work 
exists. 

IV. A. CALCULATION OF THE TANGENTIAL RADIATION PRESSURE VECTOR 
The calculation of the tangential radiation pressure vector, denoted by pr TANG , involves 
the computation of n£ r ) and (tf^- nj r ) , with 

TT er = ; Tf* = f‘<v;v + ‘> avA.ia) 


"Her 1 p° < V? V.* > 


TT+r = f'(Vv t 7 (IVA - lb> 


This obviously means that v r ''°, v 9 i o , and v^ i o must be calculated in regions interior and exterior 
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to the drop interface. This has been done in the previous section. The time averages are done 
in the same manner as Marston (1980), and are taken over an interval of 2 x/oj c . The acoustic 
carrier wave frequency is denoted by co c . It can be shown that, in the case in which the acoustic 
wave was modulated, the resulting oscillation will be at the drop frequency (with a phase factor 
also present). In the case of the unmodulated standing acoustic wave, the averaging eliminates 
any time dependence. 

From the results of the previous section it is clear that the acoustic velocity components 
are expressed as a series in "1". The r dependence is seen to involve spherical Bessel and 
Hankel functions. Also, the "inner" region stretched variables appear in exponential functions. 
The 9 dependence occurs through the expressions involving Legendre and associated Legendre 
polynomials. These functional dependencies involve infinite series. 

In particular, in nondimensional form 


(pF T "'“’) e =- (KVr'Vg > - <.V r 'v e *> (rv A. 2 a) 

Note that the velocity components are also nondimensional quantities. To lowest order (in e) 
this can be written as 
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(TVA.2b) 
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Recall that < v f0 > = < \f 6 > at the interface, and thus 
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(IVA.2d) 

However, fit - \ = 0 at the interface (based upon periodicity and continuity of pressure, 

see Marston, 1980). Therefore, 


( p?™> ) 9 



(TVA.3) 


In a similar fashion, it is possible to calculate (pr TANO ) (j ,. However, the calculation of the 
tangential radiation pressure vector will only involve the e 0 component, as axisymmetry has 
been assumed. Even so, the manipulations will involve multiplication of the aforementioned 

velocities, each represented by an infinite series. Moreover, an integration over 0 must be 

performed. 

If the arguments of the Bessel functions are small, simplifications will result. Only a few 
terms of each series then need be retained. Of course, the arguments are composed of factors 
which involve physical property data. Enough terms must be retained so as to insure that the 
solution is adequately represented. 

From inspection of equation (HA. 4c), it is clear that for any projection in which 
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the YL does not have m = 0 the integral will be zero. An evaluation of such integrals will 

involve products of expressions comprised of Legendre polynomials, with ultimately to be 
rewritten in terms of Legendre polynomials. 


46 



IV.B. SPHERICAL HARMONICS/COMPUTATIONAL CONCERNS 


In order to calculate (p? TANG ) e , it is necessary to calculate integrals such as 

1 z au $. y** { <*****>) 

(IVB.l) 

a .at . -v, 

where v ro ‘'°, v eo '’°are evaluated at r=R (f=|=0), the drop radius, and the remaining functional 
dependence is that of 0. 

The primary interest will be in the lowest order shape oscillation, and / will be taken 
to be equal to two. The resulting integrals which must be evaluated will have ^-dependence of 

the form 




(IVB.2) 


c, is a constant. 

As many / and l' values must be used as those used to represent the Bessel functions. 
Moreover, it is possible that some of the resulting integrals will be zero. 

Another form in which the tangential radiation pressure forcing could appear is in that 

given in equation (IIA. 4c). This occurs if the stress balance equations are manipulated through 

such operations as taking the surface divergence. In this case, the tangential derivatives of the 

product of, say, <v ro ‘ \£> with respect to 6 must be taken. This will be reflected in that the 
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type of terms in the integrands of the resulting integrals will become more complicated than that 
in (IVB.2). 

It is clear that this is a manipulation involving intensive calculation which is best done 
using computational resources. To that end, software development was done; the strategy is to 
make use of representations of Legendre polynomials and efficient computational schemes which 
have appeared in the literature (Press et. al., 1986). 

In particular, the integrands may be represented by Chebyshev polynomials. Also, the 
integration of a Chebyshev series can be done without quadrature per se, as relationships 
between the coefficients of the integrand and of the would be integral are known. This is 
discussed in more detail in the MSAE thesis of Ferguson. 

Of course, the Bessel functions themselves have to have been evaluated at whatever the 
constant value of the argument is. It is important to include enough terms in the series to satisfy 
that the solution is adequately represented; that is, the truncation of the terms will have to be 
done in order to obtain numerical results. 
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IV.C. EXAMPLE CALCULATION 


This section will consider the forcing in the boundary/interface conditions which will 
occur if the description of the acoustic standing wave field includes viscous terms. 

Note carefully that this section will not determine the forced hydrodynamic field. That 
is to be done in Section V and VI for the cases in which the acoustic standing wave field (with 
viscous effects included) is unmodulated and modulated, respectively. 

In a sense, this is an intermediate approach to the problem of acoustic forcing of a liquid 
drop. The work of Marston considered the acoustic field to be inviscid and the hydrodynamic 
field to be strictly represented by the natural oscillating drop field determined by Miller and 
Scriven (1968). 

The primary work of this project is to determine the viscous acoustic field and 
hydrodynamic field which is a solution to the forced governing equations. This hydrodynamic 
field will be a modification of that known to Miller and Scriven. [These are given in Sections 
III, V, and VI]. 

In this section, the acoustic field will include viscous terms. Therefore, the 
boundary/interface conditions will include tangential stresses. However, the hydrodynamic field 
will remain that of the naturally oscillating funforced') drop . 

The focus will be on the determination of the deformation to the drop due to the 
(modulated) acoustic field; referred to as the "static part" by Marston (1980). The problem will 
be done in dimensional form as to most easily compare with Marston (1980). 

It is emphasized that very few experimental measurements on a variety of drop/host 
systems have made. Therefore, calculations will be done for the case of a 1 mm drop of p- 
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xylene in water. This was the system consider by Marston (1980) and Marston and Apfel 
(1980), so some comparisons between calculations will be possible. Also, the physical property 
data is taken from the aforementioned sources. 

The other physical system with extent experimental results is a silicone oil drop in water 
(Marston and Goosby, 1985). 

The effect of parameter variation has been considered numerically (Ferguson, see MSAE 

thesis). 

The radial component of the hydrodynamic velocity field for this "static" problem is 
given (in dimensional form) by 


1 * Ci ovc.ia) 

and 

^ = J 

A* (TVC.lb) 

The Y ta (0,0) represent spherical harmonics. However, since the problem is axisymmetric, 
d 

— = 0 , and the value of m is strictly zero. The superscripts "i” and "o" indicate the 

d<p 

regions interior and exterior to the drop. 

Boundary/interface conditions are applied at r = R, with R the undisturbed drop radius. 
Note that only the radial component of the hydrodynamic field is given. In the original work 
of Miller and Scriven (1968) essentially a poloidal/toroidal field decomposition method was used 
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in the theoretical investigation. Use of surface divergence and surface curl operators on the 
boundary/interface conditions resulted in a set of relationships which involved on u/, u, 0 , and 
the radial vorticity component. (They allowed for a <f> variation.) Moreover, the conditions 
involving the vorticity decouple form those involving the velocity. (The resulting conditions 
involving vorticity pertain to the existence of shear waves.) 

In this example, only the radial (hydrodynamic) velocity components are necessary. The 
boundary/interface conditions of continuity of radial velocity across the interface and the 
kinematic condition result in 

= - C x R* (IVC.2a) 

(IVC.2b) 

as well as a consistency check. 

The continuity of the tangential component(s) of the hydrodynamic velocity is addressed 
using the surface divergence operator, with • &' - at r = R. 

Also, V n - V - V - i r — . This yields (after making use of relations is (IVC.2a-2b)) 

8r 




= 0 


(IVC.3) 


The surface divergence of the tangential stress is given by the following equation. Note that it 
is forced by the surface divergence of pr TANG , discussed in Section IVA. The constituents of 
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pr T ANG acous tic field quantities. It is 


M . 1 | aCI-rX^R 1 ' 9 R*' 1 ] 

- /j° ( -11U*D b° R' l ' H 

t 

+ / V . =1 TAWi, \ p 1 ^ * H 

\ ’'Surface j° r / K 

(TVC.4a) 

Henceforth, abbreviate V SURFACE by V s . The only component of pr TANG will be (pr TANG ) e , as there 
is no <£ dependence. Moreover, m = 0. The brackets " < > " indicate time averaging (over 
the acoustic period). 

It is understood that a summation appears in Equation (IVC.4a). Utilizing the 
orthogonility properties of spherical harmonic functions as well as Eqns. (IVC.2a-2b) yields, 
after manipulations 
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The exponential represents the phase; with the imposed phase and 17 R l the response. Let 
^dip* = Vm ' Vk- Also, via orthogonality properties, the integral which appears in the numerator 
is known to have a value of one. Therefore, Equation (IVC.4b) can be rewritten as 

(ajtm R H c i + 

* *-t( Cj> £(V S * 


(IVC.4c) 
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The final condition is that of the normal stress balance. It is given by 

<3 JJL° ^Ur “ 4* y? 

Jr 

= y a+lKt '° **.«••■> *£«(»>♦) £ n<l 

+ < . ) t ilLln 

(FVC.5a) 

As before, a summation is implied. After substitutions, manipulations and use of orthogonality 
properties, this is given by 

- V* R +' A’ 4 R l l L 

* £ (£+z)(2-0 K* lW *o 

R* 

= ( n> 

(IYC.5b) 
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The "static" deformation is given by . The surface tension/curvature term is third on the 

l.h.s., with ^indicating the surface tension. Let t? dif = - tj r ; with the phase of the 

forcing and tj R the phase of the response. 

The system of equations can be written in matrix form: 
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This is a forced system in three unknowns. The expression for k^ 0 can be obtained easily via 
Cramer’s rule. It is 


A 

Hi 




[ <r(i+i)(JH)/R*] 


( CC s * e c < v " mL > *•**) & ■' w 
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(IVC.7) 


Clearly, the tangential radiation stress will contribute to the deformation. Moreover, the 

incorporation of viscous effects into the acoustic field will modify the value of < pr radiaJ > . 

Marston calculated k^ for a 1 mm p-xylene drop in water, assuming an inviscid 
acoustic field. He took the acoustic carrier wave frequency to be 217.5 kHz. The drop’s center 

was taken to be located at the velocity modal plane of the incident carrier wave. The incident 

pressure wave had amplitude of 10 5 dyn/cm 2 . 

His calculation of pr radiJj was based on an approximate formula as opposed to expansion 

of the Bessel function series, later improved (Marston et. al., 1981). His improved formula 

calculated <pr 20 radial > as ".12 7 dyn/cm 2 . Note that 1 = 2. Deformation was found to be on the 
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0 


order of 1 micron. 

The calculation of < pr radU1 > in this work utilized the full series solution for the inviscid 
functions of y,, p, ; the velocity and pressure fields of the acoustic wave [v t ', v t °, v t \ v 9 °, p', p 0 ]; 
also projected onto the 1 =2 mode. Of course the series had to be truncated. Four terms were 
kept in the series (from j 0 to j 3 , h 0 (1) to h 3 (1) , with j, h (,) representing spherical Bessel functions). 

Physical property data was taken from Marston (1980). In this work, the incident wave 
field amplitude was expressed in terms of the acoustic velocity potential. Therefore, the value 
used by Marston for the acoustic pressure had to be converted to the corresponding value for 
the acoustic velocity. This was done. 

It was found that <pf 20 r * d “ , > equaled -. 12 dynes/cm 2 . No viscous effects were taken into 
account in the calculation of the radial pressure forcing in order that its value be compared to 
that obtained by Marston. 

In the calculation of the tangential forcing effect it was found that 

7V777 — ( S f** iL t^ r 

1 5 °] 

for 1 = 2 was equal to -0.05 dynes/cm 2 . This then modifies the amount of deformation - which 
still remains quite small relative to the drop size. 
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The above number is with set to zero. That is, the phase of the hydrodynamic field 
response was assumed equal to that of the imposed driver - which is not likely. Thus, the 
magnitude of this deformation represents a maximum . 

It is in Section V that the static deformation will be re-addressed. In that section, it will 
be shown that the newly found form of the hydrodynamic field will result in modifications 
(potentially very significant) to the aforementioned result. 

Also, in the work of this project, expansions for the hydrodynamic field were done in 
terms of Legendre polynomials, as opposed to the spherical harmonics used by Miller and 
Scriven (1968) and by Marston (1980). 

The work of Marston et. al. , (198 1) determined < pr 2 o* d “ 1 > to be -. 127 dynes/cm 2 using 
an approximate scheme. The calculation of the effect of the tangential acoustic radiation 
pressure forcing could not be determined by Marston ’s fomulation. 

Further discussion of details on the calculations of this section can be found in Ferguson 
(see Reference section). 
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V. HYDRODYNAMIC FIELD: FORCED BY AN UNMODULATED ACOUSTIC WAVE 
FIELD 


In this section, the hydrodynamic field which exists (at second order in the expansion 
parameter 8 ) as a result of the unmodulated acoustic standing wave field is investigated. In this 
case, the steady state hydrodynamic field is simply the streaming field. Moreover, the radiation 
pressure vector is comprised of static contributions only. The static radiation pressure results 
in drop deformation. This is the situation which would occur if the acoustic levitation system 
were to be used solely to position the sample rather than to induce oscillations as well. It is 
noted, however, that if the modulated standing acoustic wave field was the acoustic field, static 
deformations would exist. That is, in the case of the (induced) oscillating fluid droplet, 
deformation of the droplet will exist. Of course, if the deformation is very small relative to the 
drop dimension, it need not be taken into account in the interface/boundary conditions. 

It is remarked that equations at this order could be constructed for which the time 
dependence is exp (T 2 i t). However, they are not of interest. 

The hydrodynamic field is considered to be viscous and incompressible. It must be 
determined both interior and exterior to the drop. The generation of the second order (in 8 ) 
system of equations will be presented explicitly in regions both interior and exterior to the drop. 

VA. HYDRODYNAMIC FIELD EXTERIOR TO THE DROP: EQUATIONS AND 
SOLUTIONS 

The nonlinear system of governing equations given in Section DC (by equations IIC. 1 - 
IIC. 3) is nondimensionalized via relationships given in Section UIB (see HIB.l). 

This results in the nondimensional system 
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+ (T + t»)^. VfV U°) 

(VA.lb) 

and relationship between pressure and density which holds at order 5. Note that the "o" 

superscript refers to the region exterior to the drop. The velocity field is given by U°, the 

pressure and density by p° and p° respectively. The Reynolds-type number is denoted by Re AC , 

with the subscript denoting that quantities relevant to the acoustic field are utilized. (This has 

been referred to previously as Re). Thus, Re AC - c°(c°/co AC ]/v° . Again, the subscript on the 
frequency, « AC , indicates that it is the acoustic frequency. 

Let the dependent field variables be expanded in a series in the expansion parameter 5, 

(Recall 8 = w drop /oj ac , with a>u RO p the (natural) frequency of drop oscillation). That is, 


U° = sy,° -V ul 


(VA.2a) 
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P° = X * S P»° + 


(VA.2b) 
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(VA.2c) 


Subscripts of "1" indicate acoustic field variables, ”2" indicates a hydrodynamic field 
variable. At order 6, the governing equations of the acoustic field are recovered. These can 
be found in Section mB; see Equations (IIlB.2a-IIIB.2c). 

The hydrodynamic field equations occur at order <5 2 . Recall that the hydrodyanmic field 
is incompressible. 

Thus, 
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(VA.3a) 
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It is clear that the term (— +t"|— V /V • u°\ which occurs in the conservation of momentum 

U ) Re Ac 1 ^ 1 

equation must be zero. The quantities v, 0 and p° are known at this order, and act as forcing 
(or source) terms in the conservation of momentum equation (VA.3b). 

Interest is in the steady state streaming field. The time average of equations (VA.3a-3b) 
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over a period (2n /<o^ c ) must be taken. Before doing this, it is noted (using Equation MB. 2a) 
that 


n 


F,° !** * v,' n° 

n 'it 

f - v, # (v vH 
n 


Therefore 






^ e AC 


|r(f°v,°) + v'(vyr) 

o-t 


+ v, 0 -^\/,° 


Taking the time average yields 


(VA.4) 


(VA.5) 


(VA.6a) 

V* U° ~ O 
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( 


)l’\ (v- v,”) + y,°- VvT 

+ COMPLEX CoAJTDdATC 


(VA.6b) 


The overbar indicates complex conjugation. At this step, it is now understood that y°, pi , and 

v,° are expressions that are independent of time . It remains to solve for (and pi ), using 
system (VA.6a-6b). 

Several facets of this problem will be discussed before proceeding to a solution. The first 

aspect involves a re-nondimensionalization of the variables, and hence the system of equations. 

The second aspect will address the nature of the forcing (or source) terms. In particular, it is 

recalled from Section III that the acoustic field itself (i.e. y,°, p{ ) has been expressed in terms 
of a composite solution. The ramifications of the composite solution for this problem will be 

explored. 

Re-nondimensionalization Scheme 

Heretofore, the nondimensionalizations employed have utilized acoustic field reference 
quantities, such as w Ac and c 0 °, the speed of sound in the outer medium. The solutions which 
have be found are, of course, nondimensional. 

It is found to be more convenient at this stage to re-nondimensionalize with respect to 
reference quantities important in the hydrodynamic problem. 
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Rewriting (VA.6b) in terms of dimensional quantities yields 


1 c ) 


- _L_ JL. 7^ 


Re 


“I 


AC 



(Co/u?At} C*0 





+ V. + c 


0 


(VA.7) 


The" * " indicates dimensional quantities. 

The drop reference dimension is taken to be d, the reference velocity field is taken to 

be (dUfrop) , and the reference pressure field is taken to be (p°j (du DRO pf . These are used 
to re-nondimensionalize Equation (VA.7). 

A key point is the fact that since the acoustic frequency was chosen so 
that (c°/u> Ac } would be on the order of the drop size; the drop reference dimension d can be 

equated to (c®/wj . 

After substitution and manipulation, the re-nondimensionalized hydrodynamic field 
governing equations are found to be 


U* = 0 


(VA.8a) 
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= (-0 


V, 0 * vv,° + y,°) 

+ Complex comju&aTE 


(VA.8b) 


Note that Re 


HYDR 


d(d<D DR Qp) 
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Discussion of Forcing Terms 


The nondimensional forcing terms seen in equation (VA.8b) (or in (VA.6b)) are 
composed of known acoustic field quantities. Recall that a composite solution exists for v,°. 

A 25 

For example, v 19 ° = v 9o ° (r,0) + v 9o ° (f,0). Since the time dependence has been eliminated, this 

is a variation on results of Section III. As f -* w , v So ° -* o. There is then an acoustic "sublayer 

region" (of dimension SQRT ) outside of which any terms involving v So °, v rl °, or their 

respective complex conjugates are zero. 

Therefore, outside of this acoustic sublayer region, the only terms which contribute as 
forcing terms are of the form 

■+ Complex c<*itu&at£ 

However, expressions of this form can be re-written in terms of the gradient of a scalar function. 
The ramification is that outside of the acoustic sublayer, the forcing function, written in terms 
of the gradient of a scalar function, can be viewed as strictly a modification to the pressure field. 
Therefore, outside of the "acoustic sublayer" region, the curl of Equation (VA.8b) is 

- L V ! u” = O 

(VA.9a) 

Let x be the vorticity. Since the velocity field is two-dimensional, only 

the component of exists. Equation (VA.9a) can be written in terms of the vorticity as 
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Q 


(VA.9b) 


± ^U> z ° = 

In essence, outside of the "acoustic sublayer" region, the forcing of the hydrodynamic field is 
zero, and there are no sources of vorticity. 

It is inside the "acoustic sublayer" region only that the hydrodynamic field equations are 
forced. Solutions to this forced vector equation will be found in terms of f, the stretched 
independent variable introduced in Section HI (and 9, of course). These solutions will decay 
as f -* oo , and will represent corrections to the unforced problem in the region of the drop/host 
medium interface. 

Solution in the "Outer" Region Exterior to the Drop 

The governing equation is given by (VA.9a). It is reasonable to work this problem in 
terms of a Stokes' stream function. Let 



L 

r*$£Te 



(VA.IOa) 


O _ - t 

u ae ~ rs^e dr 

Substitution into Equation (VA.9b) yields 


(VA.IOb) 
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^ Jr* r* 3e U*& a& /J 
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L V r 1 to V** 99/ j 1 


(VA.lla) 


The solution to this equation is simply Stokes flow. It is given by 


0O 

t^Cr.e) - + F°r' i )^i^e ^(o>i&3 


(VA.llb) 


je»i 


Solution in the 11 Acoustic Sublayer" Region Exterior to the Drop 

In this region, the independent variable r is stretched, and is represented by 

r a 8 + eS 

Recall that c = 1/SQRT (Re Ac ). However, the governing system of equations has been re- 
nondimensionalized with respect to reference quantities important in the hydrodynamic problem. 

The question arises: What relationship should e have to Re^R? Recall that 5 = 
o>D R op/ ^ac* Then 




(VA.12a) 


and 


« = /s’ 


(VA.12b) 


It is remarked that other nondimensional parameters can be developed which will relate this 
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problem to that studied by Riley. Define 


R s - X 2 R Eac = - SRe H ^ 


(VA.12c) 


In this work, it is clear that Re Ac > 1 and e < 1. 

In order to select the order of magnitude of Re^DR* a relationship between the orders of 
magnitude of 5 and e must be stated. Moreover, selection of the order of magnitude of Re^DR 
will then determine the order of magnitude of R,. 

In this work, the following will be taken: 


= 00 > (VA.12d) 

Then it follows that R, = o(5). Therefore, the order of 5 must be equivalent to that of e 2 . (Of 
course, Re Ac may then be re-expressed as being of order (£').) 

A formal analogy may be made with the work of Riley. This work considered the flow 
field resulting from a solid body in oscillatory translational motion. Such a problem is quite 
different from that investigated in this project. However, there are mathematical similarities. 
In his work, the parameters Rs^ey, Cr^ey, and M,^ occur, with Rs^^ = Briley Briley- 
M riley corresponds to SQRT (Re Ac ). Criley corresponds to 5. R pn pv and Rs pn P Y correspond to 
Rchydro ^ Rs, respectively. For the case of M > 1, R pn FY of order 1, and Rs RnPY •« 1, 
Riley finds that the flow outside of a shear layer region is Stokes-like. Note that outside of the 
acoustic sublayer in this problem , the flow was found to be Stokes-flow. 
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Since e - 


_/L_ 

WTO* 


, the relationship between r and f can be re-expressed as 


r = R + yr $ 

(VA.13) 

with Rchydr of order one. It now is necessary to rewrite the governing system of equations in 
terms of the stretched variable f. Note that in this acoustic sublayer region, the forcing terms 
on the right hand side of equation (VA.8b) will contribute . These forcing terms are comprised 
of functions which involve the acoustic velocity field. Recall that this field had a representation 
in the acoustic sublayer region which involved functions of f as well as the outer region field 
(representing the inviscid acoustic field) re-expressed in terms of a Taylor series expansion. 

As was the case in the outer region, it is convenient to work in terms of a stream 
function. After taking the curl of equation (VA.8b) and expanding in terms of the acoustic 
sublayer region variable, one obtains 
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and 
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(VA.14) 
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(VA.15) 
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Inspection of the lading order fcrtns i n Equation (VA. .5) teveals that 



(VA.16) 

This simplifies to 


fit 
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*3 0a 


R SUaQ \ i 
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+ ex. 


Therefore, it is dear that is 0 f order f - 8 | in 

Then Y£, and are zero and 


the acoustic sublayer. 



It is necessary to solve for Y n . Recaii that 


(VA.17) 
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*'=l 


Let (C,0) - M) sin0 p} (cos0) , and substitute this into Equation (VA. 17). Continue the 

manipulations by multiplying through by Py (cos0) , and integrating over (o,ir). This 
eliminates the theta dependence, and yields a forced equation for f(f), which is 





with £ — / £- 


(VA.18) 



and let 
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It is found that 
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Therefore, the solution for the stream function in the acoustic sublayer region is given by 


■ X 3 (toifr) 
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* ^ ~ 0+0 s ^ 

+ (t) KaJ eyp f - (h±L$ ) 

j n * 


(VA.20) 
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VB. HYDRODYNAMIC FIELD INTERIOR TO THE DROP: EQUATIONS AND 

SOLUTIONS 

The nondimensionalized nonlinear system of governing equations is given by 



(VB.Ia) 


t w « * -V ♦ 

+ ("t — ) (VB.lb) 

and the relationship between pressure and density which holds at order 5. The "i" superscript 

refers to the region interior to the drop. The velocity field is given by I2‘ , the pressure and 
density by p' and p\ respectively. The Reynolds-type number is denoted 

b y * 

Let the dependent field variables be expanded in the expansion 
parameter h(-u DRO[ Ju A as 


u* = 

Sv,* 

+ 

(VB.2a) 

f -- 

Po 

+ + f*pT 

(VB.2b) 

p‘ = 

fl 

4 Sp ; + S l p i 
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d* Jx 

(VB.2c) 
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The quantities having a subscript "1" are acoustic field variable, "2" indicates hydrodynamic 
field variables. At order 5, the governing equations of the acoustic field interior to the drop are 
recovered. (These can be found in Section IIIC). 

The hydrodynamic field equations occur at order 5 2 - Recall that the hydrodynamic field 
is incompressible, and thus 

V 4 - O (VB.3a) 




- £* 






(VB.3b) 


It is clear that V |V • 


which appears in the conservation of momentum equation 


must be zero. The quantities v,' and p\ are known at this order, and act as forcing terms to 
the hydrodynamic field. 


Interest is in the steady-state streaming field. The time average of equations (VB.3a-3b) 
over a period must be taken. Before doing this, it is noted that 


lit: ^ 

it 
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* *JL 



(VB.4) 
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Therefore 



Taking the time-average yields 

^ • U 2 - 0 


(VB.6a) 



ot 





+ 




+ coHPLey coMjufcrt'rc 

(VB.6b) 

The overbar denotes complex conjugation. AT THIS STEP, IT IS UNDERSTOOD 
THAT pj ,and v,' represent expressions that are independent of time . 

It remains to solve for {and p£j , using system (VB.6a-6b). The methodology will 
proceed in the same manner as that of Section VA. 


Re-nondimensionalization Scheme 

As was the case in Section V A, it is of interest to re-nondimensionalize the hydrodynamic 
field equations with respect to reference quantities of relevance. 
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The reference length will be that of the drop dimension, d. The reference velocity is to 

be given by dco DROP . Finally, the reference pressure Field is denoted by p °Jd^ DROP f ■ 

A key point is that the reference length, related to the drop dimension, must also be the 

reference length for the case in which the nondimensionalization scheme is carried out with 
respect to acoustic-field-relevant reference quantities. These are tied together because the 
acoustic frequency (*> Ae ) itself was chosen so that (c^/co^) would be on the order of the drop 

dimension. Thus, d-c' c 9 Jv> Ac . 

The re-nondimensionalization of the governing equations yields 

V* U* =0 (VB.7a) 


Tip* -v -_£L Vui 

= V* • ^V,‘ ■+ V* fa y/) ) + C6WTU4MTT 


(VB.7b) 


Note that RzhydrT DROf ^l v o) • 


Discussion of Forcing Terms 

The nondimensional forcing terms which appear in Equation (VB.7b) are composed of 
acoustic field variables, which are known at this order ( 8 2 ). Recall that a composite solution was 

A • R2 . 

found for Vj*. Recall, for example, that {y : ') e = v So ‘ (r ,6) + v 5o ' (£,0). Since the time 

a ■ 

dependence has been eliminated, this is a variation on the results of Section in. As £ -» oo v 9o '4 

0 . There is then an acoustic "sublayer region" (of dimension SQRT( (v^/g> A c ) ) outside of 
which any terms involving v 9o \ v rl ', or their respective complex conjugates are zero. 

Therefore, outside of this acoustic sublayer region, the only terms which contribute as 
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forcing terms are of the form 


( ■* 1 
+ COHPIEX COMTU^ftTE - 1 


(Recall that is the acoustic velocity field potential for the inviscid acoustic field in the "outer 
region", still interior to the drop). However, expressions of this form can be rewritten in terms 

of the gradient of a scalar function. The ramification is that, outside of the acoustic sublayer, 

but still interior to the drop, the forcing function, written in terms of the gradient of a scalar 

function, can be viewed as strictly a modification to the pressure field. 

Therefore, outside of the "acoustic sublayer" region, the curl of Equation (VB.7b) is 


- « V X 1 O 

(VB.8a) 

Let (-V x u^j be the vorticity field. Since the velocity field is two-dimensional, only 
the e + component of exists. Equation (VB.8a) can be written in terms of the vorticity as 


06 

folupR 




(VB.8b) 


In essence, outside of the "acoustic sublayer" region, the forcing of the hydrodynamic field is 
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zero, and there are no sources of vorticity. 

It is inside the "acoustic sublayer" region only that the hydrodynamic field equations are 
forced. Solutions to this forced vector equation will be found in terms of |, the stretched 
independent variable, introduced in Section III (and 0, of course). These solutions will decay 
as f -* o° , and will represent corrections to the unforced problem in the "acoustic sublayer" 
region (located near the drop/host interface). 

Solution in the "Outer" Region Interior to the Drop 

The governing equation is given by (VB.8). It is reasonable to work this problem in 
terms of a Stokes’ stream function. Let 
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a© 
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(VB.9a) 
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Substitution into Equation (VB.8b) yields 


0< 

* -v 
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_ *r x 
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= 0 


The solution for \p' is simply that of Stokes' flow. It is given by 


(VB.9b) 


(VB.lOa) 
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(VB.lOb) 


oO 

^ A (r,e)- ^(A l p it3 + eV +> ) (case) 

Solution in the "Acoustic Sublayer" Region Interior to the Drop 

In this region, the independent variable r is stretched, and is given by 

r = R - el 

Recall that e-^SQRT^Re^yi . The question arises: what relationship does/should e have to 

^ e HYDR? 


Recall that 

b-u> DROF lu>Ac • Then it is seen that 



“ £ Re AC ~ $ li z 

(VB.lla) 

or, alternatively, 

€ * 

SQRT ( $ /RCh^or ) 

(VB.llb) 


It is noted that other nondimensional parameters can be developed which will relate the 
problem under consideration to the work of Riley involving a solid body in periodic motion 
(Riley, 1967). 

Clearly, the problem under consideration differs in a number of important aspects. The 
center of the drop is motionless; no solid boundaries exist in the problem. Most importantly, 
the hydrodynamic field (fluid motion) results from the existence of the acoustic standing wave 
field. The analogy with the work of Riley is that which has been presented in Section VA. 
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Define 


Rc = S X Rc AC " = S Re 


HijDR 


(VB.llc) 


In this work, it is clear that Re Ac > 1 and e <t 1. 

In order to select the order of magnitude of ReHY DR , a relationship between 5 and e must 
be established. Moreover, selection of the order of magnitude of Re^,^ will then determine 
the order magnitude of R,. 

In this work, the following will hold 


Re 


»yon 


5 0 (|) 


from which it follows that R, = o (5). 

Re-expressing e in terms of Re^DR, it is seen that 


(VB.lld) 


6 , YT_ 

Hyp** 


so that 



(VB.12) 


It is now necessary to rewrite the system of governing equations in terms of the stretched 
variable, £. Note that in the "acoustic sublayer" region, the forcing terms which appear on the 
right hand side of (VB.7b) will contribute . These forcing terms are comprised of functions 
which involve the acoustic velocity field interior to the drop. Recall that this field had a 
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representation in the "acoustic sublayer" region which involved functions of £ (representing the 
viscous correction) as well as the outer region field (still interior to the drop) representing the 
inviscid acoustic field re-expressed in terms of a Taylor series expansion. 

It is convenient to work in terms of a stream function. After taking the curl of Equation 
(VB.7b), and expanding in terms of the "acoustic sublayer" region variable £, one obtains 
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+ Complex CONJUGATE 



Inspection of the leading order terms in Equation (VB.14) reveals that 





-v Complex cowtu&ate 


(VB.15) 


This simplifies to 
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It is clear that 


Tj is of order (6/7?^^^^) in the "acoustic sublayer", such that 



(RS*» 

<* 



Vf. 




VM e , 


4 



(VB.16) 


86 



It is necessary to solve for x ¥ l 22 (i,6) . Recall that 



i I ^ . 

Let (5,0) - g(Z) sin0 P, (cos0) , and substitute this and the expression for v 6o ' into 

Equation (VB.16). Continue by multiplying Equation (VB.16) through by /^(cosO) , and 
integrating over (o, r). This eliminates the theta dependence, and renders a forced ordinary 

differential equation for g(£) which is 


- R'kswi f I'j R 1-iHijj) exf 


[%) Rii - ' ( ia i) 


with X 


4 
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and let 
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It is found that 
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(VB.18) 
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Therefore, the solution for the stream function in the "acoustic sublayer" region is given by 
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VC. BOUNDARY/INTERFACE CONDITIONS 

The boundary/interface conditions are those appropriate for the hydrodynamic field. 
Moreover, the boundary/interface conditions, applied at the drop/host interface, should include 
only linearized terms, as the hydrodynamic problem solved at this order is linear. In general, 
the conditions are: the radial component of the velocity is continuous (across the interface), the 
tangential velocity component is continuous, the tangential stress balance holds, and the change 
in normal force across the interface is balanced by the surface tension/curvature term. Finally, 
the kinematic condition holds at the interface. 

The velocity and pressure fields which enter into these conditions refer to those of the 
hydrodynamic field. (That is, these are given by pi , p 2 of Sections VA and VB. 

Of course, there are surface forces which also enter into the boundary/interface 

conditions. These are ((p? radm! )) and (ft>r TANG )l) , which have been presented in Sections II 
and IV. 

The velocity and/or pressure fields which describe the hydrodynamic field will include 
the modifications/corrections which arise as a result of the forcing by the acoustic field 
variables. Simply stated, these modifications are those contributions which arose in the "acoustic 
sublayer". Were the forcing of the hydrodynamic field governing equations to be neglected , the 
resulting hydrodynamic flow would be strictly that of Stokes flow. In the work presented in 
Sections VA and VB, it is seen that Stokes flow arises in the "outer" regions (i.e. , not in the 
acoustic sublayer), both exterior and interior to the drop. 

In the work of Marston (1980), the acoustic field was taken to be strictly inviscid. 
Then (pr TANG )^ did not contribute, and the hydrodynamic field itself was described entirely by 
Stokes flow. 

It is the intention of this section to develop the boundary/interface conditions which arise 
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in the case in which the acoustic field incorporates viscous effects. In doing this, it will be 
necessary to look at the velocity fields obtained in the previous two sections in more detail. 

As a first step, the equilibrium interface will be defined. This will be followed by a 
presentation of the general form of the interface/boundary conditions. Finally, the velocity field 
(after some manipulations) and the pressure field determined in Sections VA and VB will be 
utilized in the conditions, and specific equations obtained. Of particular interest is the resulting 
deformation (from sphericity) of the drop. 

Equilibrium Interface 

Let the equilibrium interface, Fe, be defined by 

CO 

Fe = C - R — (tose) - O (VC.l) 

t-o 

with the third term (on the r.h.s.) representing the deformation. 

Boundarv/Interface Conditions: General Form 

The general form of the boundary/interface conditions will be presented. In keeping with 
earlier work, there will have been nondimensionalized. 

They are: 

(Continuity of Radial Velocity) 



<xi r^R, 


(VC. 2a) 


(Continuity of Tangential Velocity) 
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out r-R ( S 


(Kinematic Condition) 


r= R = 

(Tangential Force Balance) 


"Ur = O a A 

V/ = O at 



(VC.2b) 


(VC .2c) 


So,<(pr w ) e ) 


ReHjjDR |(|3 Vr* J g Vj, - Vr 0 |sNe’ * CoMPliX CowTO&ftTE )j 

at r-R,^-r*o 


(VC.2d) 


Recall that and p-p^/p® . Note that and refer to the hydrodynamic field. 

The forcing due to the (0) tangential component of the radiation pressure vector is on the right 

hand side. It’s form is that discussed in Section IV. That is, for example 


a _ 
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(Normal Force Balance) 


“ P' -v 'hVf 

9 'I* 


-p° -t - G»^*A 

i <p 


(VC.2e) 

with G - (o^p" ; and o° is the surface tension associated with the interface 

between the liquid drop and the host medium. The unit vector a is the outward pointing 

normal to the interface. Specifically, a-VFe/|VFe| . 

Equations (VC.2a-2e) represent the general form of the boundary/interface conditions 
which must be satisfied. Note that forcing terms, namely, the i r and e 0 components of the 
radiation pressure vector, appear in the above system of equations. This is a nonhomogeneous 

set of linear equations. The unknowns are the coefficients of the velocity and pressure fields. 

Once these are determined, the hydrodynamic flow field which only exists as a result of the 

acoustic forcing will be known. 

Velocity and Pressure Fields 

Recall that the velocity fields in both the "acoustic sublayer" region and the "outer" 

region exterior and interior to the drop have been found (in general form) in Sections VA and 

VB. However, the stream function formulation was employed. This is not convenient, as the 

boundary interface conditions require iff and iff . Also, it is necessary to determine 
explicitly the pressure field. 

Recall that 
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oC 


t 




~** x 4. F° r'*) s^e -^(eose) 

bo + b, 5 -v b.> S 5, + b 3 
+ HKKip exp (-i^Os)l 

+ m ( Kaf ) exp (-( i - 0 ^ 


5wv6 


£ (cos©) 
(VC.3) 


However, as f-* oo, the contribution that is in the "acoustic sublayer" region must decay. This 
implies that b°- b°- b 2 - b°* 0 . 

It is a straight forward matter to determine u° 2 and uj^ . These are given by 


v/ 

S . = 




1.(0 


cp 


(B°r~* •vF°f~' g ' 3 ') 1 'i 


+ / _S \i 

R<?HUCft / R 


\i V" 7 ( (-i')(n»p6*.p(-(Ji01) 'W 9 _(sa,8R‘' 

) i‘M ^ 

+ (n(K»I)e4p(-ikLls') 
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or 



ut =. JMI+O (B°r'* + F° r'* 8 * 1 ) 

<3 i = i 

^(iXKaDe^p^^s) J 


i-\ 


d V HYDR 3 

Similarly, (recalling that — =>-* in the acoustic sublayer), 


dr 6 ac 


^0* 


- L 


rs\^9- 


00 


•'Jtl 




^ ( ’Cw^bV* + iF (tose) 
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' ( J_Y 1 V 4 

/ V|., R 


Pp (coS&) 


(VC.5) 


l + <0 (t?) ( Ki s“)^p(--^l'S )| 

Note that the contribution to from the "acoustic sublayer" region is at o^SQRT^b/Ref^^j 
and that to u‘ 2 at o^I^hydr} • Recah that Re^R is order one. 

The form of and also proceeds from tj 4 • Recall that 


tyV - ^ ( A“ r 8 * 3 + E x f *♦') s^e P^’Ctose) 

i--\ 

1d„ 4 \ \ 4 I 4 t>3 ^ 

4 (0(Ki£ ) e >P (j-JMau (1-aH) 
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(VC.6) 
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Since the contribution that is strictly in the acoustic sublayer must decay 

as S "* °°» b* 0 - b( - ~ a o • 

Then 

^r A = XL f0(|+0 (A^r^ 1 + E-V' 1 ) (cos© - ) 

ju» 



(VC. 7) 


and, recalling that in the acoustic sublayer region, — 

5r 


-A , it is found that 

5 as 


U 




- Y1 (d+v>A ; r i+l + (stotV' 1 j 


(tose) 




4 l 


Re, 




1 £ jJ vf ( 0(K '« W (- /p^ tw>0 


P j(«0») 


The specific forms of k^, u^, u r2 , u° 2 have been determined, up to a constant. The 


unknown constants which remain are F°, F®, .4*, E l . Of course, these depend upon T. 


96 



(Recall that K1 f and Klf are known from previous sections, and also depend upon "1".) 

It remains to explicitly determine p 2 and pi . From the governing equations for the 
hydrodynamic field presented in Sections VA and VB, it is clear that the contribution 

to pf from the "acoustic sublayer" region cannot be greater than order {fb IRe^^j at the very 
most. 

The determination of pf in the Stokes layer region can utilize the results of Miller and 
Scriven (1968). That is, although they solved for the hydrodynamic field of an unforced, 

naturally oscillating drop, the pressure for Stokes flow can be recovered by setting their 

oscillating frequency to zero. Then the contribution to pressure which depends on the radial 

coordinate r is proportional to, in the Stokes* flow region exterior to the drop, 


-L I r iL + 

9r \ 


< 2 . 

ar 


- 0 




This is found to be zero. A similar development can be done for the pressure in the Stokes 
region interior to the drop. 

It is then found that, to lowest order, the contribution to the normal stress balance of the 
hydrodynamic field will involve terms which are the radial derivative of the radial component 
of the hydrodynamic field. 

It is remarked that this result was also the case in the work of Marston (1980). 


Boundary-Interface Conditions: Specific Equations 


Substitution of the forms for iff, iff, p l f presented in the prior subsection into 
Equations (VC.2a-2e) will be done. 
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Lowest order terms only will be kept. It is recalled that the hydrodynamic field itself 
occurs at order 8 2 in the overall perturbation scheme. 

(Continuity of radial velocity) 

A i R , ' tl b'r" 1 + + 



(Continuity of tangential velocity) 

-(I«U'S W -»iOE i ♦ 0 

(Kinematic condition) 

B° + F° 

t £* ■* OjMfe^p*) *0 

(Tangential stress) 


We - + ±_ltU r c> ) 

Jn R R ^ 



n 



(VC.10) 


(VC.ll) 


= <(r™%> 


(VC. 12a) 


Notice that presence of the forcing term due to the tangential component of the (time-averaged) 
radiation pressure vector, but only lowest order terms will be kept. The expansion 
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for Uq (l/ 9 0 ) in terms of £ (f) can be written, utilizing a Taylor series expansion, as 


V, 


0 


(r,en 


+ /IT 




and 


+ "jit. u *. 


u, 


Ug (r, e -) t Or 

5 1 / 

+ /p <(*,») ♦ ••• 


(VC.12b) 


(VC. 12c) 


Substitution of (VC. 12b) and use of ujf- u*f 


yields (to lowest order in 


V Re 


HYDR 


), 
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Gathering terms yields 

<* A 1 (-anCi+x) k 1 <n A (-aW-on + rt R *' 1 

oe 

^ im* 4 | Ifioiel 

hi 





’‘''i-m 3 - o + f°R' J 'Y^)rt+^p;c^) 


= <(r 1 A “")s> 

Multiply through by (sin0 • P^j and integrate over ( o , t ) to obtain 

°( | A'fc* (-aiM + a^ -v £*$•*'* 

+ } * JL(k\« * wap 

« * 

- «p™%> ^6 

The right hand side of Equation (VC. 12f) can be written as 


(k>sO 


(VC.12e) 


(VC.12f) 
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Sfff+O 



using Equation (IVA.3). Note that 
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(VC.13) 


Also recall that 
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Therefore, the right hand side of Equation (VC.12f) can be further expanded as 


(3-ttO . 
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However, certain integrals are equivalent to pre-defined quantities. Therefore, the r.h.s. of 

Equation (VC.2f) can be further expressed using the definitions of and the 

continuity of the radial component of velocity of the acoustic filed across the interface. This 

yields 

(R.H.S.) of (Eq. VC. 120 




-V 




) 


Therefore, 

(-arxtta) R l A 1 + (-a’lff 2 -l)E i R 4 ' x ^ 

" t F° R '*' 3 


(VC. 14) 


\ fllj V [iXa/ - X)/ ~ 13; 


(VC. 15) 


Recall fhst 

K\; - C'lO 51 ^ 

Kl; “ t ^ 5 ( - r) Ren^Dtt 
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Therefore, the tangential stress balance in the hydrodynamic field is given as 


* [M 1 (-u)Un) + E 1 

- \_ B ° ' S '*' 1 + F ,, R ' i ‘ 3 (- jOUm)I 






\ 




(VC.16) 


Clearly, the contribution to the tangential stress balance which arose due to the modified flow 
field (different from that of Stokes’ flow) cancels the component of the forcing (radiation 
pressure) vector. This is for the very restricted subcase in which the response to the forcing has 
zero phase lag with respect to the forcing. This results in 


oM('<U)U+i) A 1 R 


* l 


-of'? 1 - 1 ] 


- 4 F 0 *' 1 ' 5 (-aOdtil) <VC ' 17 

Clearly, use of the kinematic condition could serve to simplify this expression. 

Thus, for Equations (VC. 9), (VC. 10), VC. 11), and (VC. 17) represent the specific 
boundary/interface conditions. It remains to construct the specific form of the normal force 
balance condition. 
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(Normal Force Balance) 



From the previous discussion 



Ac ( w; =r a_ 


- - <p R,m,a 


nr 


(VC. 18) 


To lowest order (in the hydrodynamic field balance) this becomes 
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(VC. 19a) 


The theta dependence can be eliminated after multiplication through via P 2 (cos0)sin0 and 
integration over (o.ir) . This yields 
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(VC. 19b) 
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Let the right hand side of (VC. 19b) be renamed as 


( (toss) -She 

n / n . . „ <2- Jn 


JUiM) 


(VC.20) 


The boundary/interface conditions applied at the order of the hydrodynamic field and 

forced by the radiation pressure vector components at the surface {r-R, C -5-o) are listed below 
for convenience. They are 


A c S i+ ' = B°R' A Vf'E’ 1 * 1 


(VC.21a) 


■(jin) A 1 R* 4 ' -(1* 0 E V ' ^ 


(VC.21b) 
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(VC.21c) 


A C ft 1 *’ + R = O 


(VC.21d) 
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and 
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(VC.21f) 


^(C+0 

The unknown coefficients are: *,, A‘, B°, E 1 , F°. (The two kinematic conditions are not 
necessary - that is, only one should be used. Either (VC. 21c) or (VC.21d) can be utilized). 

It is possible to eliminate E‘ and F°. This yields a system of nonhomogeneous linear 
algebraic equations for the unknowns A‘, and B°. These are 


R* A 1 ' + - R- 1 -' B° = 0 


(VC. 22a) 
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Final Comments 


It is clear that 



Ko = 


( I’wu,,, 


R 


(VC.23) 


K. is the coefficient of the deformation, and is directly proportional to the radial component of 
the radiation pressure vector. 

The tangential stress has been included in the calculations. First, the viscous acoustic 
field was determined. This allowed the tangential (e e ) component of the radiation pressure 
vector to be calculated. Secondly, the modification of the hydrodynamic field due to source 

terms in the governing equation was determined. The source terms are caused directly by the 

inclusion of viscous effects on the acoustic field. Solution of the governing equations shows 
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there to be enhanced flow in the "acoustic sublayer" region, particularly of the u© component. 

I ft ) 


The flow correction due to this enhanced flow is o 


for Uq, which is the largest such 

\ i^HYDR J 

Stokes* flow is the flow in the "outer" regions, 


component in the acoustic sublayer region, 
both interior and exterior to the drop. 

In the application of the boundary/interface conditions, only lowest order contributions 
(at the level of the hydrodynamic field) are included. Of course, in the tangential stress balance 
condition, this results in the inclusion of velocity gradient terms (of Uq) from the acoustic 
sublayer region. 

It was found that these contributions to the stress balance cancel the forcing due to the 
tangential component of the radiation pressure vector for the special subcase in which there is 
no phase lag of the hydrodynamic field with respect to the initial forcing. 


It is clear that incorporation of viscous effects into the acoustic field and investigation of 
the forced hydrodynamic field shows that the velocity field has a level of complexity which is 
not apparent if only the inviscid acoustic field is considered. It is only through doing the 
"viscous" problem that the structure of the velocity field in the "acoustic sublayer" regions - 
which border the drop/host interface can be elucidated. The hydrodynamic (streaming) field 
due to an unmodulated acoustic standing wave field has been considered. In the case of a 
modulated standing wave field, static (streaming) contributions to the hydrodynamic field would 
exist as well as the oscillating contributions. However, there would be a phase lag between 
forcing a response, in general, so that the deformation would be due not only to the radial 
component of the radiation pressure vector, but to the tangential component as well. 
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VI. HYDRODYNAMIC FIELD: FORCED BY A MODULATED ACOUSTIC STANDING 
WAVE FIELD 


This section details the structure of the hydrodynamic field which exists (at second order 
in the expansion parameter 6 ) strictly as a result of the modulated acoustic standing wave 

field. In this case, the resulting hydrodynamic field is oscillatory in time. The drop itself is 
undergoing (small amplitude) shape oscillations due to the forcing at the natural frequency of 
the drop. 

Let (a AC -(j3 , -(ji ,, +<j3 DRO p , with a" the acoustic frequency of a second acoustic wave. 

Clearly, , the drop oscillation frequency. The frequencies of the acoustic waves 

are orders of magnitude larger than that of the natural frequency of the drop. The time- 
averaging of the governing equations over the period [ 1 will result in a set of time- 

dependent equations. However, the frequency of the resulting time dependence is at the "beat" 
frequency, which is that of drop’s natural frequency of oscillation. (Refer to Miller and Scriven, 
1968, for a discussion of natural (free) drop oscillations.) 

The hydrodynamic field is considered to be viscous and incompressible It must be 
determined both interior and exterior to the drop. The generation of the governing equations 
for both regions will be presented explicitly. 

It is remarked that the methodology of this Section is quite simialr to that of Section V. 
Of course, the time dependence introduces an additional level of complexity. However, some 
of the discussions, such as the need for re-nondimensionalization, apply in both cases. Where 
possible, it is the intention to avoid reduplication providing clarity is not compromised. 
Moreover, an expanded discussion as to the motivation behind several of the methological steps 
will be included. 


ill 




VIA. HYDRODYNAMIC FIELD EXTERIOR TO THE DROP: EQUATIONS AND 
SOLUTIONS 

The nonlinear system of governing equations given previously in Section IIC (see Eqs. 
IIC. 1-3) is nondimensionalized through utilization of the relationships given in Section IIIB (see 
IHB.l). Acoustic field reference quantities were used in the nondimensionalization scheme. 
The resulting sytem is 



4 V*. + 




= O 


(VIA. la) 



4 p °u 0 -yu 0 = ~sjp° + j^v 2 V° 


+ 


v[V- V°) 

v ^ / Reac. 


(VIA. lb) 


Also, the relationship between pressure and density holds at o ( 6 ) . The "o" superscript refers 
to the region exterior to the drop. The velocity, pressure and density fields are given 

by u°, p°, and p° , respectively. Re AC is a Reynolds-type number, and is given as 
c o{co/^ A cj/^l .with (j) ac the acoustic frequency. 

Let the dependent variables be expanded in a series in the expansion parameter 6 . 

(Note that &-u> DR0P / u> AC .with <* DROP representing the natural frequency of oscillation of the 
drop.) 

That is, 
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(VIA. 2a) 


V* S\j,° 4 £*uj 

r ■- c w ♦ *v 

f* = 1 + # + s V 


(VIA.2b) 


(VIA. 2c) 


Acoustic field variables are denoted by the "1" subscript (y^, p°, p°) , and hydrodynamic 

quantities by "2" (u%, p 2 °) . At order 5, the governing equations of the acoustic field are 
recovered (see Section IIIB, Equations IIIB.2a-2c). 

It is at order 5 2 that the hydrodynamic field occurs. It owes its existence to the 
(modulated) acoustic standing wave field. The hydrodynamic field is taken to be incompressible. 

Resulting governing equations are then 


^ uj = o 


(VIA. 3a) 


- i \x{ 

D-t 
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(VIA. 3b) 



Clearly, the term ( — +r°) — - — V /V • u°) which formally appears at this order is zero due 

to the incompressibility of the hydrodynamic field. 

The quantities y° and p° are known at this order. They act as forcing (or source) 
terms in the conservation of momentum Equation (VIA. 3b). 

Now; the acoustic field represented by Y° (and, of course, p° ) is a modulated 
standing wave field. For example, 

0 1 ' ' 1 " 

:i = + (VIA.4) 



with (the scattered) velocity potential for the acoustic wave. The refers to the acoustic 


wave of frequency g> / (-(d ac ) and """ refers to the wave of frequency 


The time average of Equations (VIA3a-3b) over a period of 



must be taken. This 


will result in having a time-dependent source term (i.e., the right hand side), of a form to 
be exhibited explicitly. A solution for ( and p 2 °) will be sought with this same time- 
dependence. Taking the time average of terms on the right hand side yields 


<RHS> 


j>* 1 1 + <*<u>s(sjb +y-y) 

( tv*) | Si + a tos (Si *'»| w 

+ ( Y,° • Vv,° )( A* -y)) 

* ($r • ^ a -y) 
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(u 




Me) 


(VIA.5) 
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In this expression only, the caret " A " indicates a time independent general quantity. Both 
r\" and V indicate phases. Note that there is a steady state contribution to the forcing as well 

as the cos ( 5 fc+V /_ V ) time dependent forcing. The steady state forcing was addressed via 
consideration of the unmodulated acoustic standing wave field, in Section V. It will not be 

considered in this section. Rather, the focus will be solely on the response of the hydrodynamic 

field to the oscillatory forcing. 

The governing equations, after being time-averaged, are 

^ ~ ° (VIA. 6a) 




‘Pa 


- L V 1 




AC 


Vl 1 MV,° -v 


1 g* ( t ^\° ^ 

-V COMPLEX COUT06.ATE 


(VIA. 6b) 


co +> 1 " 'l' ) 

(S - oft op /oJac ) 

Note that Equations (VIA.6a-6b) are nondimensional, and have been nondimensionalized 
with respect to acoustic field reference quatitiies. The periodic temporal forcing contains "5t"; 
with 6 <? 1. However, at this stage, the hydrodynamic field is of interest. A re- 
nondimensionalization scheme will be performed following the methodology discussed in Section 
V. 
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Re-nondimensionalization 


From a physical point of view, it is the hydrodynamic field which is the focal point. It 
is the hydrodynamic field which exists only as a result of the forcing. This is clearly shown in 
Equations (VIA.6a-6b). However, in these equations all quantities have been 
nondimensionalized with respect to reference quantities which are not relevant to an oscillating 
drop. 

The re-nondimensionalization acts as a renormalization, with the result that the oscillating 
drop problem is viewed from the standpoint of the hydrodynamic field as opposed to that of the 
acoustic field. 

Rewrite (VIA. 6b) in terms of dimensional quantities as follows 

/— +• / 1 Vfp/') 

\fi ] ii 
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c) C-o] V ♦ £ (ivD 

<f^ 

-V COMPLEX COMT06.MTE 


• COS ( S CJac ^ - vj ’ ) 


(VIA .7) 


116 



The " " indicate dimensional quantities. Renondimensionalize as follows: 


$ 






Mx 5 x 


<Llo 




OftOp = kk 


(AOozc*) X P = V 5 


(VIA. 8) 


The reference length is that of the drop dimension, d. This, in fact, provides the coupling to 

Q ° 

the acoustic field; d — . The natural frequency of oscillation of the drop provides the time 


b> 


AC 


scale. 


Utilizing the re-nondimensionalization scheme yields, after manipulations 


V -%1 =0 


(VIA.9a) 
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+ f t ° ( 

+ COMPLEX CeWTUSAfC 


k>S ( ± 


(VIA.9b) 

with Re^ jj- (d 2 w d /v£) . The quantities are now nondimensional with respect to 

hydrodynamic field quantities. Moreover, Si? and 9° represent the hydrodynamic/drop 
oscillation field quantities. The system represented by (VIA.9a-9b) is linear. Of course, the 
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right hand side of (VIA. 9b) contains nonlinear forcing terms, but these are known quantities. 
It proves to be more convenient to work in terms of 


Uz ~ ° (VIA.IOa) 

dt ^e^DR. 

= C-ai ^ ^ (VU - 10b) 

Let 

* hrl ' U°(r,^,9) (jfc “'‘l 1 ) ) 

■^4 “ f a ° (VIA. 11) 

Substitution of (VIA. 11) into system (VIA.lOa-lOb) yields 

^ ~ 0 (VIA. 12a) 
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f ”uvn +c .c.) (vu ub) 


It is the system (VIA. 12) which must be solved for (and p 2 ° ). 


Discussion of Forcing Terms and Resulting Decomposition 

The forcing terms on the right hand side are of the form 
* Vj£ + pi (i^) ♦ c.c.j , and represent real quantities. In particular, 

V? * v°(r,B) + 

A 

with £° referring to the velocity field in the outer region (represented by V<J>° ), exterior to 

the drop and with ( ( , 0 ) referring to the correction to the acoustic velocity field in the 
"acoustic sublayer" region. As f -* oo, quantities which depend on f will go to zero. 

If the curl of Equation (VIA. 12b) is taken, the quantities on the right hand side which 

contain only terms representing the inviscid outer flow will not contribute. That is, if the 

vorticity equation is constructed, it will be unforced in the outer region. 

With this in mind, decompose as follows 


Ua = & (f,e) + u° ( $,e) 

fsT = p/M) + 

(VIA. 13) 
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Solution in the "Outer" Region Exterior to the Droi 

Taking the curl of Equation (VIA. 12b) - and the restricting attention to the "outer" region 
exterior to the drop (i.e. no remaining C dependence), one obtains 

- _t - o 

3 t few,,* ~ (VIA. 14) 

with the vorticity in the outer region exterior to the drop. Notice that there are no sources 

of vorticity; in this region the hydrodynamic field is unforced. However, this is precisely the 
problem that was investigated by Miller and Scriven (1968), although their analysis was 
dimensional and the notation differs somewhat from the present exposition. 

The solutions can then be written for both the velocity and pressure fields in this region. 


A ° , 
Ur,tf,e) 
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(VIA. 15a) 
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(VIA. 15c) 
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with a 2 - (-i R e HYDR) - 

Solution in the "Acoustic Sublayer 11 Region Exterior to the Drop 

In this region, the independent variable r is stretched, and the independent (radial) 
variable C is related to r as 

r - R + rs 


Recall that e- - — . However, the system of governing equations has been re- 

nondimensionalized with respect to hydrodynamic field reference quantities. 

1 

A relationship between e ( - ) and Re H ^ must be stated. 

Re ac 

Recall that 8 - w MOP /w AC . Then 

= £(R*ac) = (VIA. 16a) 


£ ~ / V Ke H r 

a (VIA. 16b) 
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One can also define a third "Reynolds number-like" parameter, R g , such that 


R S ° = SRe H;jW (VIA. 16c) 

In this work, it is clear that Re AC » l and e < l , 6 < 1 . If the order of magnitude 
relationship between 5 and e is specified, an order of magnitude for R, will follow. 

The remarks concerning the relationship of this forced drop problem to the work of Riley 

(1967), which involved a solid body in periodic motion which were presented in 

Section V carry over to this Section. Therefore, they will not be repeated here. 

Select Re HYDR - 0(1) . Since e-yfS/ y JRe HWR , the relationship between r and 5 can 
be re-expressed as 



(VIA. 17) 


with Rchydr of order one. It is now necessary to rewrite the governing system of equations in 
the variable f. Note that in this region, the forcing terms on the right hand side of (VIA. 12b) 
will contribute . The forcing terms themselves are comprised of functions which involve the 
acoustic velocity field. 

The expansion will be done in terms of the primitive variables as opposed to a stream 

function. Since the conservation of momentum equation is a vector equation, the 
<§ r and S g component equations will be expanded in the acoustic sublayer region. Writing 

the S r and e 0 component equations yields: 
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The right hand side forcing terms in (VIA.18a-b) (written terms of f) are known from 
knowledge of y ° (and p£ ), which are the acoustic field quantities. 

Therefore, 

( e r cons, momentum - expanded) 
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The Y° and p° quantities represent the acoustic field. They are expanded in the 
acoustic sublayer region; and contain both terms which decay to zero out of the sublayer region 

(i.e. as f -* oo ) and terms which represent the inviscid acoustic field quantities re-expressed in 

a Taylor series expansion in terms of f in the sublayer region. The subscript "1" 
on y° and pi is kept in order to distinguish clearly these acoustic field quantities from the 
hydrodynamic field quantities, which are indicated by a "2" subscript. It should be clear that 

there is not any time dependence present in Equations (VIA.18a-18b) through (VIA.20a-20b). 

For more details concerning the acoustic field quantities, see Section III. 

The expanded form of the e e conservation of momentum equation will be exhibited. 
The above remarks also apply to this equation. 

( e 0 cons, momentum - expanded) 
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Note that in both (VIA.20a-20b), the terms on the right hand side involve products, as least one 

factor which involves only the viscous correction to the acoustic field. Also, it is necessary to 
expand the conservation of mass equation in order to relate u° 2 and Qe° . It is 

(cons, mass-expanded) 



(VIA.21) 

It is clear from Equations (VIA.20a-20b) that a balance can be achieved between the 
forcing terms and the hydrodynamic field terms if is selected to be of order ■ - ^ . 

\^ e HYDR 

Physically, the hydrodynamic field is forced by the acoustic field. The objective of the 
mathematics is to have the lowest order hydrodynamic field terms balanced by the lowest order 
forcing terms. 
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Inspection of (VIA. 20b) shows that this is the case if the highest derivative of the 
tangential component of the hydrodynamic field balances the radial derivative (normal to the 
interface) of the tangential component of the acoustic field, convected by the radial component 


of the inviscid acoustic velocity field. 


Moreover, it must be that is zero, and that the 


balance occurs with £3^ (£ , 0) . 


This has ramifications for the u£f. component. From Equation (VIA. 21), it is clear 

that % is the lowest order contribution to %■ 

From (VIA. 20b), 
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Then, using (VIA.26) in (VIA. 25), and multiplying through by 
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, and integrating 
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yields, 
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which yields 
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The solution to (VIA. 27(1, can be written as 
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VIB. HYDRODYNAMIC FIELD INTERIOR TO THE DROP: EQUATIONS AND 
SOLUTIONS 

The nonlinear system of governing equations (see Section IIC) is nondimensionalized with 
respect to acoustic field reference quantities. The resulting system, governing flow interior to 
the drop is: 



•V 
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= o 


(VIB. la) 


s n 


(VIB.lb) 



with a-p^p" and t'-p^g/p® . The "i" superscript refers to the region interior to the drop. 
(Also, the relationship between pressure and density holds at order 5 .) The velocity, pressure, 
and density fields are given by IT, p l , and p‘ , respectively. Re AC is a Reynolds-type 

number, and equals c^{c^/(a A( ^lv o 0 , with u AC the acoustic wave frequency. 
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Let the dependent variables be expanded in a series in the expansion parameter, 6 . 


(Note that ^“^dhopI^ac ^drop indicating the natural frequency of oscillation of the 

drop.) That is 

<$v + S X Ua 

(VIB.2a) 






(VIB.2b) 


f - ^ + 


(vm.2c) 


Acoustic variables are denoted by the "1" subscript, (v,', p,', p,'), and the hydrodynamic 
quantities by "2" p£j . At order 5 , the governing equations of the acoustic field are 
recovered (see Section III). 


It is at order 5 2 that the hydrodynamic field occurs. It owes its existence to the 
(modulated) acoustic standing wave field. 

The hydrodynamic field is taken to be incompressible. Resulting governing equations 
are then 


V • u‘j =o 


(VIB.3a) 
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Clearly, the term ~+x ! 


Re 


— V /V • u^j which formally appears at this order is zero due to the 

AC 


incompressibility of the hydrodynamic field. 


The quantities v,‘ and are known at this order. They act as forcing or source terms 
in the conservation of momentum Equation (VIB.3b). 

The acoustic field represented by v,‘ (and p(, p\ ) is a modulated acoustic standing wave 
field. For example, 



» A 9 

4 > + 



(VII!. 4) 


with the velocity potential for the acoustic standing wave. The refers to the acoustic 

wave of frequency , and """ refers to the wave of frequency m" . 

The time average of Equations (VIB.3a-3b) over a period of 2 tt/w j4C must be taken. 
This will result in having a time-dependent source term (i.e., the right hand side), of a form to 

be exhibited explicitly. A solution for (and P[ ) will be sought with this same time 
dependence. 

Taking the time-average of terms on the right hand side yields 
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-V- p ( vf • vv; ) ( £+ <2 6 0 s(£±*>)"-f )) 

i p>(yF'V)>i ) (a + aeosfa+y-y)) 

S z k>p«op /^ac. (VIB.5) 

In this expression only, the caret indicates a time independent, general quantity . 

Both r\" and t/ indicate phases. Note that there is a steady state contribution to the forcing 

as well as the cos(6/ + r| // - r] / ) time-varying dependent forcing. The steady state forcing was 
addressed via consideration of the unmodulated acoustic standing wave field (see Section V). 

It will not be considered in this section. Rather, the focus will be solely on the response of the 

hydrodynamic field to the oscillatory forcing. 

The governing equations, after being time-averaged, yield 
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1 ' ) (VBB.6b) 


Note that Equations (VIB.6a-6b) are nondimensional, and have been nondimensionalized with 
respect to acoustic field reference quantities. The periodic temporal forcing contains "5t", with 
5 « 1 . However, at this stage, the hydrodynamic field is of interest. A re- 
nondimensionalization scheme will be performed following the methodology discussed in Section 
V. 


Re-nondimensionalization 

The hydrodynamic field is the focal point. It is this flow which exists as a result of the 
acoustic forcing. This is clearly shown in Equations (VIB.6a-6b). However, in those equations, 
all quantities have been nondimensionalized with respect to reference quantities which are not 
relevant to an oscillating drop. In particular, the "5t" in the periodic time dependence 
indicates a very slow scale - yet the drop oscillation time should define the time scale of the 
hydrodynamic field. 

The re-nondimensionalization acts as a re-normalization, with the result that the 
oscillating drop problem is viewed from the standpoint of the hydrodynamic field as opposed to 
the acoustic field. 

Rewrite (VTB.6b) in terms of dimensional quantities as follows 
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The " A " indicates dimensional quantities. Re-nondimensionalize as follows: 
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(VIB.8) 

The reference length is that of the drop dimension. This, in fact, provides the coupling to the 
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acoustic field : d - c°l<ji AC . The natural frequency of the drop oscillation provides the time 
scale. 

Utilizing the re-nondimensionalization scheme yields, after manipulations 


V- 'U* X 0 


(VIB.9a) 
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• <U)S( 1 + >(/' (VIB.9b) 

with RChydr = (d 2 oj d /V 0 °). The quantities are now nondimensional with respect to acoustic field 

quantities. Moreover, Ti^and Pj, represent the hydrodynamic/drop oscillation field quantities. 
The system represented by (VIB.9a-9b) is linear. Of course, the right hand side of 

(VIB.9b) contains nonlinear forcing terms, but these are known quantities. It proves to be more 

convenient to work in terms of 

~ 0 (VIB.lOa) 
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Discussion of Forcing Terms and Resulting Decomposition 


The forcing terms on the right hand side are of the form jpy| + p + c.c.j , 

and represent real quantities. In particular, the field variables yj and p, are functions 
of (r,0) and £ . For example 

v; = V* (r,6 1 

with v,' referring to the velocity field in the outer region interior to the drop and with j,' 

(£,0) referring to the correction to the acoustic velocity field in the "acoustic sublayer" region. 
As £ -» oo , quantities which depend on £ will decay to zero. 

If the curl of Equation (VIB. 12b) is taken the quantities on the right hand side which 
contain only terms representing inviscid outer flow will not contribute. That is, if the 
(linearized) vorticity equation is constructed, it will be an unforced equation in the outer region. 
With this in mind, decompose as follows: 


~ ^a(r,e) * 

" Pa'lr.a) t P;u,e) 

(VIB.13) 

Solution in the "Outer" Region Interior to the Drop 

Taking the curl of Equation (VIB. 12b) - and restricting attention to the "outer" region 
exterior to the drop (i. e., no remaining £ dependence), one obtains 
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(VIB.14) 


with 4 ^ the vorticity in the outer region interior to the drop. Note ther are no sources of 
vorticity; in this region the hydrodynamic field is unforced. 

However, this is the problem investigated by Miller and Scriven (1968), albeit their 
analysis was dimensional, and differed somewhat from the present exposition. 

The solutions can then be written down for both the velocity and pressure fields in this 
region. They are: 
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Solution in the "Acoustic Sublayer 11 Region Interior to the Drop 

In this region, the independent variable r is stretched, with the independent (radial) 
variable ? related toras P = R — £ 5 • 

Recall the e = l/SQRT(Re AC ). This system of governing equations has been re- 
nondimensionalized with respect to the hydrodynamic field reference quantities. A relationship 
between e and Re^n must be stated. Recall that 



£ (*e AC ) - <£ /& 2 ’ 


(VIB.16a) 


and 


e - yrr 

One can also define a third Reynolds like number, R,, with 


(VIB.16b) 
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£ I 


(VIB.16c) 


Rs = = 

In this work, Re AC > 1, and e < 1, 5 •« 1. Now Re^DR is taken to be order one, and so R, 

< 1 . 

The remarks concerning the relationship of this forced problem to the work of Riley 
(1967), which involved a solid body undergoing periodic motion; which was presented in Section 
V carry over to this section. They will not be repeated here. 

Since ReHY DR = 0(1), e = 5/ Re^oi* is then of order SQRT (5). The relationship 
between r and £ is given by 



(VTB.17) 


with RChydr of order one. It is now necessary to rewrite the system of governing equations in 
the stretched variable £. Note that in this region, the forcing terms on the right hand side of 
(VIB.12b) will contribute . The forcing terms themselves are comprised of functions which 
involve the acoustic velocity field. 

The expansion will be done in terms of primitive variables as opposed to the stream 
function. Since the conservation of momentum equation is a vector equation, the e r and e„ 
component equations will be expanded separately in the acoustic sublayer region. Writing the 
6 r and e„ component equations yields: 
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Now, expand as 
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A similar expansion is done for pj'. All functions overbarred with a "~ n are dependent upon 
£ and 8 . The right hand side forcing terms in (VIB.18a-18b) (written in terms of f) are known 
from knowledge of v,' (and p,'), which are acoustic field quantities. 

Therefore, 

(e r cons, momentum - expanded) 
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The v,' and p,' quantities represent the acoustic field. They are expanded in the acoustic 
sublayer region; and contain both terms which decay to zero as £ -* oo and terms which 
represent the inviscid field quantities re-expressed in a Taylor serves expansion in terms of £ in 
the sublayer region. The subscript "1" on v,' and p,' is kept in order to distinguish clearly these 
acoustic field quantities from those of the hydrodynamic field, indicated by a "2" subscript. 

Note that the terms on the right hand side involve products, at least one factor of which 
involves only the viscous correction to the acoustic field. It should be clear that there is no 
time-dependence in Equations (VIB.18a-18b) through (VIB.20a-20b). For more details on the 
acoustic field quantities, see Section III. 

The expanded form of the t 6 - conservation of momentum equation is needed at this 
point. It is 

(e 9 cons, momentum - expanded) 
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Finally, it is necessary to expand the conservation of mass equation in order to relate 
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(VDB.21) 

It is clear from Equations (VIB.20a-20b) that a balance can be achieved between the 
forcing terms and the hydrodynamic field terms if is of order ( 5/ Rc^dr). 

Physically, the hydrodynamic field is forced by the acoustic field via convective 
acceleration type terms; that is, the forcing terms are quadratic in the acoustic field contribution. 
The objective of the mathematics is then to have the lowest order hydrodynamic field terms 
balanced by the lowest order forcing terms. 

Inspection of (VIB.20b) shows that this is the case if the highest derivative of the 
tangential component of the hydrodynamic field balances the radial derivative (normal to the 
interface) of the tangential velocity component of the acoustic field, convected by the radial 
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component of the inviscid acoustic velocity field. Moreover, it must be that is zero, and 
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that the balance occurs with 
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This has ramifications for the component. From Equation (VIB.21), it is clear that 


* is the lowest order contribution to u^. 
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which, when the right side is written out, is 
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(YIB.23) 
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After a substitution of (VIB.26) into (VIB.25), multiplying through by (d P^/d0)sin 0, and 
integration over (o,x), 
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The solution to (VIB.27d) is found to be 
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VIC. BOUNDARY/INTERFACE CONDITIONS 


The boundary/interface conditions are those that must be applied at the hydrodynamic 
field level. Of course, these conditions themselves will be linearized, as the hydrodynamic 
problem itself is linear. The boundary/interface conditions which must be applied are: (a) 
kinematic condition, (b) continuity of velocity field components across the drop/host medium 
interface, (c) tangential stress balance across interface, and (d) normal stress balance (including 
surface tension/curvature contribution) across the interface. 

The velocity and pressure fields which are utilized are those of u 2 ', u 2 °, p 2 i , p 2 °, given in 
Sections VIA and VIB. 

In addition, the surface forces due to the radiation pressure vector will enter into the 
boundary/interface conditions. Recall that these are the projection of the acoustic radiation stress 
tensor upon the surface of the drop. (That is, pr^ 01 ^ and pr TANG , which have been presented 
in Sections II and IV.) 

The velocity (and pressure) field which comprise the hydrodynamic field do include 
contributions which arise as the result of forcing (of the governing equations) by acoustic field 
quantities. These contributions arise in the acoustic sublayer. If their contribution to the 
hydrodynamic problem velocity field were to be neglected, the resulting hydrodynamic flow field 
would be only that of an oscillating drop (previously investigated by Miller and Scriven, 1968). 
In the work presented in Sections VIA and VIB, it has been seen that the flow in the "outer" 
regions (both interior and exterior to the drop) is precisely that of an oscillating drop. Of 
course, in the work presented in these sections, the contribution to the velocity field arising in 
the acoustic sublayer region is determined, also. 
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In the work of Marston (1980), the acoustic field was taken to be strictly inviscid . Then 
a non-zero contribution to (pr TANa )* does not exist. Moreover, the hydrodynamic governing 
equations themselves are not forced. Therefore, only (pr" dU1 ) will force the oscillating drop 
solutions in the boundary/interface conditions in the case in which the acoustic field is strictly 
inviscid. 

It is the goal of this section to develop the boundary/interface conditions which apply in 
the case in which the acoustic field incorporates viscous effects. In doing this, it will be 
necessary to look at the velocity fields obtained in the previous two sections in more detail. 

As a first step, the equilibrium interface will be defined. This will be followed by a 
presentation of the general form of the interface/boundary conditions. Finally the velocity and 
pressure fields pertaining to the forced hydrodynamic problem will be utilized in the conditions, 
and specific equations obtained. 


Equilibrium Interface 

The equilibrium interface is defined by 
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The third term (following the first equality) represents the time-dependent oscillation of the 
interface. 


Boundary/Interface Conditions: General Form 

The general form of the boundary/interface conditions will be presented. In keeping with 
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earl,er w °*, these have been 


nondimensionaJized. They are: 
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Recall that a = p 0 7p 0 ° and /S = p 0 7p 0 °. U l r and U l r refer to the time-independent factor in 
the hydrodynamic field. The forcing due to the tangential component of the radiation pressure 

vector is shown (on the r.h.s.). Its form is that discussed in Section IV. The " A " over^ r J° 

refers to the contribution due only to the inviscid acoustic field. Clearly, \f„ 0 i o , the£„ component 

of the acoustic field, includes viscous effects. That is to say, for example, 
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with G = (oJp° d 3 a> DR0P 2 ), and a° the surface tension associated with the interface between the 
liquid drop and the host medium. The unit vector fi is the outward pointing normal to the 
interface. Specifically, n = V F d | VFe | . 

Equations (VIC.2a-2e) represent the general form of the boundary /interface conditions 
which must be satisfied. Note that forcing terms, namely, the e r and e 6 components of the 
radiation pressure vector, appear in the above system of equations. This is a nonhomogeneous 
set of linear equations. The unknowns are the coefficients of the velocity and pressure fields. 
Once these are determined, the hydrodynamic flow field, which exists as a result of the acoustic 
forcing, will be known. 


Velocity and Pressure Fields 

Recall that the velocity fields in both the "acoustic sublayer" region and "outer" region 
interior and exterior to the drop have been found - up to the unknown coefficienct - in Sections 
VIA and VIB. Also, the pressure field in the "outer" region (interior and exterior to the drop) 
was determined explicitly (also up to a constant.) The pressure field in the "acoustic sublayer" 
region is of higher order in (y/5/ Re^o*) and will not contribute. 

Recall that 
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However, as f -* oo, the contribution that is in the "acoustic sublayer" region must decay. This 
implies that d 0 ° = d,° = 0. Therefore, 
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The pressure field exterior to the drop is given by 
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Note that the contribution to the pressure is essentially due to the "outer" region solution. 

The velocity and pressure fields interior to the drop must be utilized also. These are 
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The requirement that the solution which involves the "acoustic sublayer" variable, £, decay as 
£ -* oo implies that d 0 ‘ = d/ = 0. Therefore 
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Finally, the pressure field is given by 
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Boundary Interface Conditions: Specific Equations 


It is necessary to substitute the expressions for u^ 0 , u^°, p^ 0 found in Equations 
(VIC.3-VIC.8) into the general form of the boundary/interface conditions given by Equations 
(VIC.2a-2d) in order to obtain the specific system. 

Lowest order terms only will be kept. 

(Continuity of Radial Velocity) 
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(Continuity of Tangential Velocity) 
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(Kinematic Condition) 
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(Tangential Stress Balance) 

This requires some further manipulation in order to see clearly the substitution steps. 
Expanding on Equation (VIC. 2d) yields 
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Of course, after substitution, only the lowest order terms are kept. This leads to 
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But 


= 4 . 


# 4 £^( 5 ^ (S A t a;, 4f»*j 

' 

" (f Sjl 


C i ** S’ ) 


L + complex ^ujo^te 


a e * ,it Y , °) , > 



£ tt t* * £ c* a 

LM 1=1 da 


/ fiO 

-(S<S 


^ + 2 ?: 

da l*i 


at\ 

Ja / 


(VIC. 13a) 


(VIC. 13b) 


dP 2 

Therefore, multiplying by — - , and integrating over (o,x) yields, after utilizing orthogonality 
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properties, 
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However, recall that 
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Therefore, taking ^ - JL 
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The factor [21(1+1)/(21+1)] appears in each term (and can be cancelled). This yields 
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Note that the contributions due to the velocity field gradient in the acoustic sublayer, which are 
"2 jS RChydr (cl,* + c2,')" on the left hand side and "2 Re^R (cl,° + c2,°)" on the right hand 
side are exactly cancelled by the tangential component of the radiation pressure vector which 
appears as "2 j3 Re^DR (cl,* + c2,') - 2 Re,*^ (cl° + c2, 0 )" on the right hand side! 
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Therefore, the net effect is that the tangential stress resulting from the velocity gradient 
in the acoustic sublayer region variable cancels the tangential component of the radiation 
pressure vector at the interface. Both of these effects owe their existence to the inclusion of 
viscous effects in the acoustic field. 

As far as this particular boundary/interface condition is concerned, there is no change 
from the case in which the acoustic field itself was taken to be strictly inviscid. 

It must be noted that the phase of the response has been set to be equal to the phase of 
the forcing. If there were to be a phase lag, then the exact cancellation would not occur . 

Rather, let ^imposed = v" - v' and let ^response = Vr be unspecified. In the case in which 
Vr — imposed i total cancellation occurs. However, if t? r ^ Imposed (= tjoJ, then the tangential 
stress balance will be given by 
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That is, if there is a phase lag between the forcing and the response, then the tangential 


stress balance becomes a forced equation (in unknowns (a 2 °, a/, a, 0 , a 3 ')), with 
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with (pr radiai ) the time independent contribution of the radial component of the radiation 
pressure vector. 

Notice that in this section, the possibility of a phase lag between the forcing term and the 
response is acknowledged explicitly. 

After substitution, the normal force balance equation becomes 
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It is understood that there is a summation sign (over 1) in the above equation. Recall that a 2 = - 
i RChydr and fi 2 = -i ROhy DR (/3/a). 

Multiply (VIC. I9b) through by (P<£ sin 8) and integrate over (o,x). Use of orthogonality 
properties yields 
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Clearly, this can be simplified to 
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with* — -i R^hydr <tnd A -i RChydr 08/o). 

Looking more closely at pr rsdUi , it is seen to be the projection of the radial component of 
the radiation pressure vector onto the "P^" mode (or mode, if spherical harmonics required). 
It is of the form 
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A. 

with p 1,0 the adiabatic compressibility, and the subscript "1" indicating the acoustic field. Now, 
to this order of approximation; v° remains the same as if only the inviscid acoustic field were 
taken into account. It is in v te° that the viscous correction to the acoustic field will contribute. 

Therefore, although the form of the forcing term is the same as in the case in which the 
acoustic field is inviscid, it is in fact "corrected" by inclusion of the terms due to viscous 
acoustic field contributions. 
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It is remarked that the drop is taken to be spherical. The actual deformation to the drop 
was determined in Section V. Calculations indicated that for a millimeter drop, the deformation 
from sphericity is on the order of microns; and is therefore neglected as a possible modifying 
agent in the boundary/interface conditions. 

Final Comments 

It is helpful to the purpose of further discussion to reprise the final form of the 
boundary/interface conditions in one subsection. They are given by 
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Clearly, this is a forced system in five unknowns: a£ a 3 ^, a£ a*", and If the forcing terms 
are set equal to zero, the problem reduces to that of free drop oscillations addressed by Miller 
and Scriven (1968). 

In comparing, it must be noted that they did not nondimensionalize. Also, they found 
it more convenient to work in terms of surface divergence conditions. (However, their results 
can be recovered from those listed in Equations (VIC.21a-21e), provided the forcing terms are 
set equal to zero, the system re-dimensionalized, and some manipulations performed.) 

If there is no phase lag between the time periodic forcing and the time periodic response; 
then only the normal force balance equation is forced. However, even in this case the results 
do differ from those in which the acoustic field is considered strictly inviscid - that is, the 
viscous acoustic field would modify the forcing in Equation VIC.21e). 

To summarize, the inclusion of viscosity in the acoustic field produces three effects: 

1. ) in the "acoustic sublayer" region, the velocity field is enhanced 

V 

2. ) the tangential stress balance equation, taken at the drop/host interface, r = R, £ = f = 

0. (linearized case), is forced by the (projection of the) tangential component of the radiation 
pressure vector provided there is a phase difference between the time periodic forcing and the 
time periodic response. 

and 

tp 

3-) which occurs in the forcing of the normal stress balance is modified by inclusion 

of viscosity in the acoustic field. 
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vn. CONCLUSIONS 


The focus of this project has been to elucidate the nature of the flow field (both interior 
and exterior to the drop) for the case in which acoustic radiation pressure forcing due to the 
tangential radiation stress is incorporated. 

It is clear from the physics of this problem that viscous effects in the acoustic field must 
be included in order to have a non-zero tangential radiation pressure contribution. Prior 
analytical work in this area is represented by the major work of Marston (1980), and several 
other related later papers (Marston et. al., 1981, 1982). In this prior work, the acoustic field 
was assumed to be inviscid. Thus, no tangential stresses could be taken into acount due to the 
nature of the formulation. 

Physically, it is recognized that the flow field of the oscillating drop only exists if the 
acoustic standing wave field itself is modulated, with the difference between the acoustic wave 
frequencies that of the natural oscillation frequency of the drop. In the mathematical formulation 
of this situation, prior work has noted that, were a consistent formulation to be developed, there 
would be no forcing of the hydrodynamic field vorticity equation since the acoustic field was 
taken to be inviscid. Therefore, the previous analytical formulation assumed that the 
hydrodynamic flow field would be adequately represented by the known solution of the freely 
oscillating drop; and simply proceeded to force this solution by adding the forcing term the 
boundary/interface conditions. 

In addition to the primary goal of this research project which did determine the viscous 
acoustic field and elucidate the nature of the resulting hydrodynamic field forcing, certain side 
topics were investigated briefly. These include the radial forcing only of a compound drop with 
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core/shell of the same density and compressibility (presented in Appendix II) and a discussion 
of the effects of non-axisymmetric acoustic forcing. It was seen that shear waves could be 
excited were the acoustic forcing to be non-axisymmetric. 

Work done in this project has determined the viscous correction to the acoustic field. 
These results are presented in Section III. It is remarked that in order to accomplish this task, 
an acoustic boundary layer/method of composite expansion technique was used. 

The acoustic boundary layer/composite expansion approach was developed as an 
alternative to an exact anlytical method approach involving an acoustic field decomposition that 
was tried initially. The results of this first approach are listed in Appendix I. Briefly, this first 
approach was not useful because it was too difficult to implement the resulting 
boundary/interface conditions numerically - some terms proved to be extremely small relative 
to others in the equation. It is this result which pointed in the direction of a boundary 
layer/composite type expansion approach. 

Results obtained for the acoustic field with viscous effects (Section III) included show that 
the solutions for the "inner" region decay exponentially as the inner region variable goes to 
infinity. An inner region solution was found for both regions interior and exterior to the drop. 
Since the drop is liquid as opposed to a solid sphere, the perhaps more typical case of imposing 
the no-slip condition is replaced by solving a set of boundary/interface conditions. One 
interesting result of this process was that the boundary/interface conditions of (1) the radial 
component of the acoustic velocity field balance at the interface and (2) the normal force balance 
at the interface were unchanged from what they had been for the case in which the acoustic field 
is strictly inviscid. (Of course, in the work of this project, nondimensionalizations have been 
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done. This must be remembered in any comparisons with the literature.) The inclusion of 
viscous effects in the acoustic field is manifest in the boundary/interface acoustic field conditions 

of (1) tangential (<? 6 ) component of the acoustic velocity field balance at the interface, and (2) 
the tangential stress balance at the interface. It was determined that the tangential stress balance 

is due only to terms which arise via incorporation of the viscous effects. In the tangential 

\ 

velocity balance, both terms which arise due to viscous effects (i.e., the inner region 
contribution) and to the inviscid formulation (i.e., the outer region contribution) enter. It is 
through the solution of these boundary /interface conditions that the unknown coefficients are 
found, and the acoustic field characterized. 

Of the inner region dependent field quantities, only the tangential component of the 
acoustic velocity is of order one (in the inner region, it decays to zero as the inner region 
independent variable goes to infinity). All others of lower order. 

The results of Section III which give the viscous correction effects to the acoustic field, 
although interesting in their own right, are absolutely necessary in order to determine the 
hydrodynamic field. It is in Sections V and VI that the structure of the hydrodynamic field itself 
is elucidated. 

Before addressing the hydrodynamic field structure, the ramifications of incorporating 
viscosity into the acoustic field description on the tangential radiation stress tensor and on the 
radiation stress vector were explored in Section IV. The form of the tangential radiation stress 
was exhibited. 

It must be kept in mind that few examples in the literature exist on real or even supposed 
experimental systems. The work of Marston and others has utilized p-xylene drops in water and 
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silicone oil drops in water. The calculations of relevant quantities in the inviscid acoustic field 
approximation - are minimal in the existing literature. A caluclation due to Marston (1980) and 
later improved/corrected by Marston et. al., (1981) for a p-xylene drop in water exists. It is 
this result which forms the basis for a quantitative comparison. The work in Section IV has 
proceeded to calculate the forcing terms in the boundary interface conditions only . This is an 
intermediate step which may be done before actually solving for the new hydrodynamic field 
itself. 

The comparison between the calculations of this project for a one mm drop of p-xylene 
in water and those of Marston (1980, 1981) were done in Section IV. The results between the 
work of this project and Marston ’s were found to be quite close for the calculation of the radial 
component of the radiation pressure vector, which was on the order of (.1) dynes/cm 2 for a 
carrier wave amplitude (in pressure) of 10 5 dynes/cm 2 . This comparison utilized terms found 
in the inviscid acoustic field approximation only. In the calculation of the effect of the tangential 
stress upon the drop deformation, it is necessary to utilize the viscous acoustic field. It was 
found that effect is of the order of between (.01) to (0.1) dynes/cm 2 . The formulation of 
Marston (1980, 1981) could not and does noj include a viscous acoustic field, and so did not 
calculate this quantity. 

It is possible to do a limited parameter study of "theoretical" systems. Because no 
comparisons exist with experiments, the results are not being presented here, but will be found 
in the MS thesis (by Ferguson). 

The results of Section IV show that the viscous acoustic field contributes tangential 
stresses which have an effect upon drop deformation, as well as other quantities. 
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It is the actual structure of the hydrodynamic field which has been probed in Sections V 
and VI. 

In the work of Marston, the flow field is taken to be that of a freely oscillating/decaying 
viscous drop which is forced via the boundary/interface conditions. No modification of the 
hydrodynamic field itself occurs due to the forcing of the Navier-Stokes equation by functions 
quadratic in acoustic field quantities. This entirely consistent with the level of approximation 
done by Marston; but such an approach cannot incorporate tangential forcing, which is the 
primary goal of this project. 

A formal expansion scheme in the small parameter 6 , which represented 

{^DRorl^ acoustic) was utilized in obtaining the forced hydrodynamic field incompressible 
Navier-Stokes equation. 

It is noted that a re-normalization was performed at this stage. This involved a re- 
nondimensionalization with respect to hydrodynamic field reference quantities as opposed to 
acoustic field reference quantities. Following this, the following restriction was made that 

Re HYDRO (“ U DROP ^ o) is ° f 0rder 0ne - 

A detailed discussion relating this to streaming Reynolds number as well as Re AC was 
done in Sections V and VI. The relationship of this work to that of Riley was explored. 

The resulting forced equations were solved. It was seen that the forcing of the equations 
by acoustic field quantities only existed in an acoustic sublayer region; outside of this region the 
flow field is given by either Stokes flow (if the acoustic field was unmodulated) or the freely 
oscillating drop flow field of Miller and Scriven (if the acoustic field was modulated). 

This, then, represents the modification in the hydrodynamic field due to acoustic field 
forcing. Moreover, is is seen that the correction to the hydrodynamic velocity field is smaller 
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than order one in the acoustic sublayer region. Moreover, these contributions decay as the 
"acoustic sublayer" region is left. It is the tangential component of velocity that is the largest 
in this sublayer region, and therefore the tangential component of the hydrodynamic field which 
receives the most enhancement from the acoustic forcing. However, this largest enhancement 
is smaller than order one. 

The application of boundary/interface conditions at this hydrodynamic field level produces 
a very interesting result. 

The boundary/interface conditions at the hydrodynamic field level are the kinematic 
condition, the velocity component balance at the interface, and the stress balances at the 
interface. It is the normal force balance across the interface which includes the surface 
tension/curvature contribution. These quantities are all evaluated at the interface. 

There are also contributions due to the radiation pressure vector; the radial component 
of which forces the normal force balance equation. It is the e $ tangential component of the 
radiation pressure vector which forces the tangential stress balance equation. 

The radiation pressure terms have an associated (constant representing the) phase, and 
the hydrodynamic field variables have an associated phase; in general these are not equal. 

The effect of the viscous acoustic field contributes most dramatically to the tangential 
stress balance. Not only do the tangential velocity components resulting from acoustic field 
forcing in the sublayer appear in this equation, but also the tangential component of the radiation 
pressure vector serves to force this equation. Of course, the tangential velocity components 
associated with the standard oscillating drop solution also contribute. 

After the evaluation at the interface, it is seen that the velocity field modification terms 
due to the acoustic forcing in the sublayer are the same but opposite in sign from the terms 
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involved in the tangential radiation pressure forcing - up to the phase angle. If the phase of the 
response (i.e., the hydrodynamic field) is the same as the phase of the forcing, then the 
contributions would identically cancel! This would result in the effect of the acoustic forcing 
being present only through the radial component of the radiation pressure vector. (Of course, 
this radial component itself is recalculated to include the viscous contributions of the acoustic 
field, but the form would not change from that used by Marston.) 

However, it is not reasonable in general to suppose a response (i.e., the hydrodynamic 
field) which does not have a phase lag with respect to the forcing. The smaller this lag, the 
more the tangential forcing term is reduced . The final system given in Section VI may be solved 
for values of the physical parameter of interest. 

Finally, although there is enhancement of the flow field in the "acoustic sublayer" region, 
this enhancement is less than order one. Also, outside of this sublayer region, the flow is that 
given by Miller and Scriven (1968), and assumed by Marston (1980); i.e., the flow field of the 
oscillating drop. 
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APPENDIX I : 


DETAILS ON VELOCITY DECOMPOSITION APPROACH 
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This section presents in detail the determination of the viscous acoustic field by means 
of the velocity decomposition approach. 

The acoustic field which incorporates the effects of viscosity must be solved for in the 
regions interior and exterior to the drop. Quantities in the respective regions will be designated 
by superscripts i and o respectively. 

To facilitate solution of the acoustic field equations, the velocity vector is decomposed 


as 

V," = vi j) -1 + vxA"* 


(Al.la) 


v,” = v4>° + 

into irrotational and solenoidal parts. The vector potential is denoted by A' (and A°). In order 
to render that analytical calculations tractable, axisymmetry is assumed. Then A 1 = (0,0,^), 
with ^(r,0,t) a scalar function which is independent of <£ This is done in a similar fashion for 
A°. The vector potential has a structure akin to that of a toroidal field in the familiar 
poloidal/toroidal vector filed decomposition approach. However, \p does not generate A in the 
usual sense. 

It is noted that exterior to the drop the acoustic field must account for the scattered as 
well as the incident wave contributions. (This is in addition to whatever contribution exists due 
to the presence of viscosity.) Thus, <t>° = 4>£, + <f>» ie°- 

Substitution of the velocity field decomposition into the linearized governing equations 
and further manipulations yields 


VxA e 


(Al.lb) 
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Exterior region 




6a 


31’ 


(tf f I 1 *±J i< ♦ A»l 1 ) V 1 4>1, 

1 m » f.* i at I 


= 0 


(AI.2a) 


!£. = 

frt (AI.2b) 

Interior region 




1 + 



3 ^ I h] 


= o 


(AI.3a) 


(AI.3b) 

The incident wave is known strictly in terms of ^c 0 . In order to detemmine the scattered wave 
field, both A, 0 and <t>,° must be determined. This is also the case when determining the velocity 
field in the interior of the drop. 

Note the factor which appears (with i,o superscripts) in Equations AI.2a and AI.3a. This 
indicates that the effects of (shear) viscosity (ji) and bulk viscosity (/x B ) appear in the 


f 


a* V* A 
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longitudinally propagating part of the acoustic wave. These effects serve to attenuate the wave. 
Moreover, the bulk viscosity does not contribute to the solenoidal part of the velocity field. This 
is evident by considering Equations AI.2b and AI.3b. These equations are vector diffusion 
equations, with the shear viscosity responsible for the diffusion. This will be discussed further. 
The incident traveling wave can be expressed as 



% 

A lMt ^ e | 


t-o 


(AI.4) 


where w is the acoustic frequency (carrier frequency in the case of the modulated wave). This 
wave is proceeding from -oo to the drop in the direction of the polar axis. Ultimately, the 
interest is in the acoustic standing wave field. This can be expressed as the superposition of two 
(oppositely) traveling wave fields. 

Assuming a form of <t>,° (or <f>') to be an expansion in Legendre polynomials, with 
harmonic time dependence 

oo 




« 0 ''p 

^ Pje Cee»o) 4> c 




Jt*o 


(AI.5) 


the equation governing the longitudinally propagating part of the wave reduces to 


AO A 

+ (K 0 l - Jt(M ) 4-“ - O 

Jt<x r dr r* 
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with 



fa !cj ) * 


l l '«pU 




P 




(AI.6b) 


A change of superscripts will denote the equation governing 4 >\ The solution for the exterior 
region is given by 


Ct fa 9 ' 0 ^ ^en»)e' JUi 

J e*o (AI.6c) 

with h, (1) (K 0 r) indicating the outwardly propagating wave. In the interior of the drop, the 
solution is given by 





(AI.6d) 


Since axisymmetry has been assumed, there will be only one contribution to the vector 


potential equations, the e* component. This yields a scalar equation 



(AI.7a) 


governing the shear wave exterior to the drop. Replacement of the o by i superscripts, together 
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with the physical properties data, gives the equation governing the shear wave interior to the 
drop. The exterior solution is given by 


C Cr ' B '° * £ 


-leak 




(AI.7b) 


where 7 = SQRT (i w p o °/p 0 °) and is the associated Legendre polynomial. 
The interior solution is given by 


SO 

J fa ( V. r) 


1 . N -JLU>~k 


(AI.7c) 


M 


Since the standing wave field is of interest, these solutions will be superimposed with 
waves travelling in the other direction. This can be indicated by a change in the multiplicative 
constants. Let 


* (0* fat+O^ 


(AI.7d) 


with h the distance from the acoustic velocity nodal plane to the drop’s center. (This will be 
taken to be zero; then, for odd values of 1, 5, will be zero.) 
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with 


V* -- st a;* = si a° 

•- At : S t A i 

The general form of the solutions which satisfy conditions of finiteness in the respective 
regions has been determined. It remains to solve for the multiplicative constants, which can be 
done through the imposition of the appropriate conditions at the drop/host medium interface. 
The conditions on the acoustic field are that (1) the pressures balance, (2) the radial and 
tangential components of velocity are continuous (respectively), and (3) the shear stress balances 
at the interface. This yields four conditions for each 1 value, and there are four unknown 
constants existing for each 1 value. It remains to determine the pressure in order to construct 
all the necessary conditions. 

The pressure can be determined quite readily. It is a solution of 


V'° = - of* 

91 


(AI.9a) 


The pressure in the exterior and interior regions is given by 

p'Cr.s.tl ~ j <>£(&* + rar -^TST 


Xu /t 


I^Aiuc * YW ~~$r) j f* 


(AI.9b) 


p ‘foe,*)* (- 




(AI.9c) 
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For completeness, the velocity components interior and exterior to the drop are listed 
below. These are 








£. P e 
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)) it(ltr) t 


‘Ctot (Al.lOa) 




(ALlOb) 




UtSri 


^TstT^ 0 )^ “ 


+ A s * C '° r fr^ r k (Wir 0 $1 t 

JL 

X *% ^(**0 )^ (ai.ioc) 

* + i: < u; 4 (^>)^v wp ^ e '" wl 

H* &*sU'V^) 4» 6 (Al.lOd) 

+ S.£ ^‘)4*4r ( r AJVor)) fje («>s*)e ' xa * 

The conditions at the interface are then given by 

Ki i^«cKcio\ * 4 ^oitfrcR) 

5 ^ Ko (KoR^ ■* Asj 
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( *t - aJ |(<10 jtttiK) - ii ( tcR^ 




‘ A*. \ (r„o ' r° Aui+i + Alwc ^(KqR - ) 

R l R 1 R (Al.llb) 


R x 

+ 3' «*}!<, (*<Dj 

= ^rfjajOjpttaO-iCtHyo] 

+ As * t 411 *"' 5 * H?R'] Jibuti + a&RAt+iOuol 

(AI.llc) 

O^t (*f*CoVi*>) | ~ 4 ) ^jt(KiR) “ 

■t < (AK) S [. Kin) - KiV] |t(K;R') + i«<« jWK.R')} 

R x 

+ Aji (t-lljiffiO - V t Rj f+l U £ R')|| 

- ** UflWV«l J t(»*oC'/A»-Ol5 , (K.O-'«? - Jlt2tto«| 

+ ^jc l(lti)-ft?R 1 ’l -Ajt’CKiR') t ati»Wi+, (K.R'i] 

f A °* { *•"$'<** ' nmaiwifl (AI Ud . 
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APPENDIX H; 

ACOUSTIC FORCING (UNMODULATED) OF COMPOUND FLUID DROP: 
SUBCASE OF CORE/SHELL VISCOSITIES DIFFERENT 
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This work considers the effect of acoustic forcing upon a compound drop system. The 
compound drop system itself consists of a core fluid of density p 0 \ viscosity pj, surrounded by 
a shell fluid (density p 0 \ viscosity p 0 ") embedded in a host medium (with p 0 °, p 0 °). In general, 
the surface tension values at the inner and outer interfaces will differ. Only a specialized sub 
case of the general compound fluid drop system is considered. Attention is restricted to the case 
in which the core and shell fluids have the same density and compressibility, but differ in 
viscosity . 

The acoustic wave providing for the forcing of the drop is considered to be inviscid. 
That is, the viscous effects will be restricted to the hydrodynamic field. The acoustic wave then 
sees the outer interface as a boundary between two media since the acoustic wave is taken to be 
inviscid, the only contributing term in the radiation pressure vector will be the radial component. 

The case of the unmodulated acoustic wave is considered. This will result in drop 
deformation. The problem is to determine the deformation and compare it to what would exist 
for a simple drop. 

Linearized governing equations for the hydrodynamic field reflect the steady state nature 
of this problem. Manipulations on the conservation of momentum equation will yield the 
following system, expressed in the radial components of velocity, vorticity (of the hydrodynamic 
field) 

(ruj \) ii{ i.s.o] 

^*(rto r 0 * 0 ^ 
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(AH. la) 


(AD.lb) 




That is, the steady state forced hydrodynamic field problem can avail itself of the same 
decomposition as that used by Miller and Scriven in their analysis of the unforced oscillating 
drop problem. The boundary/interface conditions can be expressed strictly in terms of the 
aforementioned radial components. 

Since the only forcing component from the acoustic field is in the radial direction, the 
necessary boundary conditions will be expressed in terms of the radial velocity only. Solutions 
to the system are found to be 


^Cr,e,4))s 7T 


£,m 




(ADL2a) 




,4) - J ( M w + tit *)%,**) 




(AH.2b) 


Urtr.e.+’l - X ( Wl 1 '*'* + (Ur' 1 ) V t „ (», 

«,« 


♦ ) 


(AH.2c) 


The boundary/interface conditions which must be imposed are the kinematic condition, 
continuity of velocity components across the interface, and that the normal and tangential stress 
balances at the interface. This is at both the inner and outer interfaces. However, the additional 
(acoustic) forcing in the radial direction will enter into the normal stress balance at the outer 
interface. It is this contribution which will force the system, and result in deformation. 
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Imposition of the bounds conditions ieads to the fo«owin g systenn of options 
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With s = SQRT (R,R 0 ) and r = SQRT (R./R,). This system of equations is a set of forced 
equations in the unknown coefficients: c,, a,, b„ c„ d„ d 0 , f , f°. The deformation at the outer 
interface is indicated by f°. Note that this is a reduction from the tenth order algebraic system 
which would have been expected if the compound drop were oscillating (either free or forced). 
This is due to the repetition of relationships in two interface equations. 

This is an (8 x 8) nonhomogeneous linear system which can be solved by standard 
numerical means. Although not efficient numerically, it is clear that an application of Cramer’s 
rule would lead to a solution for f°, which is of interest. This solution would be in the form of 
a numerical value. Application of the symbolic mathematics manipulator REDUCE (Hearn, 
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1987), leads to an expression for f°, which is 



t'(hi)U-i) 




(AH.4) 


This is the result obtained for the deformation of the interface of the simple drop. Note that the 
inner surface tension does not contribute. At first, this appears to be a surprising result. 
However, upon reflection it is recalled that the core and shell region fluids differ only in their 
respective viscosities, and that the unmodulated acoustic wave responsible for the forcing is 
taken to be inviscid. 
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